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Abstract 


The Langley 14- by 22-Foot Subsonic Tunnel’s data reduction 
software utilizes six major functions to compute the acquired data. These 
functions calculate engineering units, tunnel parameters, flowmeters, jet 
exhaust measurements, balance loads/model attitudes, and model /wall 
pressures. The input( required) variables, the output) computed) 
variables, and the equations and/or subfunction) s) associated with each 
major function are discussed. 


Introduction 

In the Langley 14- by 22-Foot Subsonic Tunnel, data acquisition consists of analog, digital, 
pressure, and automated cart data. Analog and digital data is acquired with the NEFF 600, pressure data is 
acquired with the ESP 8400, and automated cart data is acquired using the Channel Access software 
package. After data acquisition, data reduction functions are utilized to compute engineering units, tunnel 
parameters, flowmeters, jet exhaust measurements, balance loads/model attitudes, and model /wall 
pressures. For consistency, the names written to an online output data file are used as much as possible in 
this document. All calculations are single precision except where noted. 
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Engineering Units’ Conversions 

Analog channel data is acquired and converted to millivolts. Next, the data is converted to 
engineering units. The engineering units’ conversions commonly used in the 14- by 22-Foot Subsonic 
Tunnel for analog data are as in, lin, poly, and type _j. 


ASIN 

* Equation Definition 

value = ( (chan#### - chan#### woz ) - bias ) / sensitivity 
if (value >1.0) then 
value =1.0 
if (value < -1 .0) then 
value = -1.0 

eu = ( dasin(value) * 57.2957795131 ) - offset 

* Setup Example 

name,FALPHA chno,3 dau,Neff 
desc, Model pitch 
units, deg 
euconv,asin 

lce(l. 233, 260.45, -0.102) 

range, 160 

filter, 1 

cmpllmt,-10 

cmpulmt,10 

dztype, no 

* Equation Example 

value = ( chan0003 - 1 .233 ) / 260.45 
if (value >1.0) then 
value =1.0 
f (value < -1.0) then 
value = -1.0 

FALPHA = ( dasin( value) * 57.2957795131 ) + 0.0102 
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LIN 


Equation Definition 

eu = ( ( chan#### - chan#### woz ) 

* ( expected power supply voltage / 

actual power supply voltage ) * sensitivity ) + offset 


Setup Example 1 

name,TA chno,l dau,Neff 

desc, tunnel ambient temperature 

units,degf 

euconv,lin 

lce(0.036,-40) 

range, 5 120 

filter, 1 

dztype, no 

Equation Example 1 

TA = ( chanOOOl * 0.036 ) - 40.0 


Setup Example 2 

name,NF738 chno,10 dau,Neff 

desc, Body Balance Normal Force 

units, lbs 

euconv,lin 

lce( 14 1.1 5448,0) 

filter, 1 

cmpllmt,-400 
cmpulmt,400 
dztype, wind 
pscorr(V738,5) 

Equation Example 2 

NF738 = ( ( chanOOlO * ( 5.0 / V738) ) - 

(chan0010 woz * ( 5.0 / V738 woz ) ) ) * 141.15448 
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POLY 


Equation Definition 

eu = constant 1 + ( (chan#### - chan#### woz ) * constant2) 

+ ( (chan#### - chan#### woz ) 2 * constant3) 
+ ( (chan#### - chan#### woz ) 3 * constant4) 
+ ( (chan#### - chan#### woz ) 4 * constant5) 
+ ( (chan#### - chan#### woz ) 5 * constant6) 
+ ( (chan#### - chan#### woz ) 6 * constant?) 
+ ( (chan#### - chan#### woz ) 7 * constant8) 

Setup Example 

name,SROLL chno,43 dau,Neff 
dcsc,IRTS Shaft Roll 
units, deg 
euconv,coef 
poly, 3 

lce(0. 1847,-0.039,-0.2865,0.3239) 

range, 160 

filter, 1 

cmpllmt,-0.5 

cmpulmt,0.5 

dztype, no 

Equation Example 

SROLL = 0.1847 + ( chan0043 * (-0.039) ) 

+ ( ( chan0043 ) 2 * (-0.2865) ) 

+ ( ( chan0043 ) 3 * 0.3239) ) 
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TYPE_J 


❖ 

❖ 

* 


❖ 


* note 

* note 


Equation Definition 

type J denotes iron-constantan 

value = ( chan#### + jcn ) 

* execute subfunction tcic(value, temp ) 
eu = temp 

Setup Example 

name,ESPTMP52 chno,65 dau,Neff 
desc,Temp of ESP52 
units,degf 
euconv,typeJ 
range, 10 
filter, 1 
dztype, no 
jcn,0 

Equation Example 

value = chan0065 

* execute subfunction tcic(value, temp ) 
ESPTMP52 = temp 


type_k denotes chromel-alumel; execute subfunction teea 
type_t denotes copper-constantan; execute subfunction tcee 
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Tunnel Parameters 


Input Variables 

ccust - customer dynamic pressure calibration constant (default = 0.0) 
dpicut - dpils/dpihs cutoff parameter (psf) 

dpihs - differential pressure of sidewall static pressure in the entrance cone 
referenced to ptot, used for high dynamic pressure (psf) 
dpils - differential pressure of sidewall static pressure in the enhance cone 
referenced to ptot, used for low dynamic pressure (psf); 
currently, value set to dpihs during engineering units’ calculations 
kprcoef - dynamic pressure probe coefficient 
pa - ambient pressure (psf) 

pitot - test section dynamic pressure as measured by pitot tube (psf) 

ptot - settling chamber total pressure (psf) 

qcode - dynamic pressure computation code; constant 

qcstn - dynamic pressure station; constant 

rjg - gas constant for air times gravitational constant 

ta - entrance cone ambient temperature (degF) 

tdew - dew point temperature (degF) 

wcode - test section configuration code for classical wall corrections; constant 

Output Variables 

dpi - measured indicated difference between total and static pressure (psf) 

dpinf - dynamic pressure uncorrected for compressibility (psf) 

mach - free stream mach number with blockage corrections 

machu - uncorrected free stream mach number 

mu - absolute viscosity of air corrected for temperature (lbs/sec) 

pstat - test section static pressure with wall corrections (psf) 

pstatu - uncorrected test section static pressure (psf) 

pv - vapor pressure calculated from dew point (psf) 

q - tunnel dynamic pressure corrected for compressibility; wall corrections (psf) 
qu - uncorrected tunnel dynamic pressure with compressibility (psf) 
rho - air density corrected with blockage corrections (slugscf) 
rhou - uncorrected air density (slugscf) 

rnft - test section reynolds number with blockage corrections (1/feet) 
rnftu - uncorrected test section reynolds number (1/feet) 
tr - ambient temperature (degR) 

vel - test section free stream velocity with wall corrections (ft/sec) 
velu - uncorrected test section free stream velocity (ft/sec) 

* note Reference 14- by 22-Foot Subsonic Tunnel document entitled “Dynamic Pressure 

* C alibration Corrections” for the current derivation ofcprime. 

* note Cprime( dynamic pressure calibration coefficient) is calculated using equations 

* determined by utilizing Check Standard Probe test results. 

* Cprime is not written to output data file 
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* default values 


dpicut = 20.0 
rjg = 1716.2290 
kprcoef = 0.998 
pitot = 0.0 

* tunnel total temperature (degrees rankine) 

tr = ta+459.67 

* vapor pressure 

pv = 12.7654*exp(((9.72334*tdew)-3 1 1 . 147)/ 

((0.555556*tdew)+223. 192)) 

* tunnel air density 

rhou = (ptot-pv*0.3789)/(rjg*tr) 

* air viscosity 

mu = 0.0002672/(tr+198.72)*(tr/518.69)**(1.5) 

* QCODE Calculations 

if (qcode > 0) then 
dpi = 0.0 
dpinf = 0.0 
if (qcode >1) then 
if (dpihs > dpicut) then 
dpi = dpihs 

else if (dpils > O.O.and.dpils <= dpicut) then 
dpi = dpils 
endif 
endif 

if (qcode = 1.0) then 

* Test Section PITOT Probe 

pitot = amaxl(pitot,0.0) 
pstatu = ptot-pitot 

else if ((qcode = 2.0) and (qcstn = 1.0)) then 

* Closed Test Section, Boundary Layer Removal System Suction is Off 


dpinf = cprime*kprcoef*dpi 
pstatu = ptot-dpinf 
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else if ((qcode = 3.0) and (qcstn = 1.0)) then 


* Closed Test Section, Boundary Layer Removal System Suction is On 

dpinf = cprime*kprcoef*dpi 
pstatu = ptot-dpinf 

else if ((qcode = 4.0) and (qcstn = 1 .0)) then 

* Open Test Section, Boundary Layer Removal System Suction is Off 

dpinf = cprime*kprcoef*dpi 
pstatu = ptot-dpinf 

else if (qcode = 5.0) then 

* Open Test Section, Boundary Layer Removal System Suction is On 

dpinf = cprime*kprcoef*dpi 
pstatu = ptot-dpinf 

else if (qcode = 6.0) then 

* Customer Assigned Calibration Constant 

dpinf = ccust*dpi 
pstatu = ptot-dpinf 

else if (qcode = 1 1.0) then 

* Walls Up, Ceiling Down, Boundary Layer Removal System Suction is Off 

* This calibration is out of date. Need new calibration for future test(s). 

dpinf = cprime*kprcoef*dpi 
pstatu = ptot-dpinf 

else if (qcode = 12.0) then 

* Walls Up, Ceiling Down, Boundary Layer Removal System Suction is On 

* This calibration is out of date. Need new calibration for future test(s). 

dpinf = cprime*kprcoef*dpi 
pstatu = ptot-dpinf 

else if (qcode = 13.0) then 

* Walls Up, Ceiling Down, Boundary Layer Removal System Suction is Off, 

* Vortex Yanes(Generators) are On 

dpinf = cprime*kprcoef*dpi 
pstatu = ptot-dpinf 
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else 


pstatu = ptot 
endif 


* note compressibility correction 

* gamma = 1.4 (gas constant); -1.0 / gamma = -0.7143 

* density 

if (pstatu > 0.0) then 
rhou = rhou*(ptot/pstatu)**(-0.7143) 
else 

rhou = 0.0 
endif 

* mach number 

if ((pstatu > 0.0) and ((ptot/pstatu) >= 1.0)) then 
machu = sqrt(((ptot/pstatu)**(2. 0/7.0)- 1.0)*5.0) 
else 

machu = 0.0 
endif 

* dynamic pressure 

if (wcode = 0.0) then 
xlamd = (ptot-pstatu)/ptot 

qu = 3.5*ptot*((1.0-xlamd)**(5.0/7.0)-1.0+xlamd) 
else 

qu = 0.7*pstatu*machu*machu 
endif 

if (qu < 0.001) then 
qu = 0.0 
endif 

* air velocity 

if (rhou > 0.0) then 
velu= sqrt(abs(2.0*qu/rhou)) 
else 

velu = 0.0 
endif 

* reynolds number 

rnftu = rhou*velu/mu 


else 


9 


* values for qcode = 0.0 

dpi = 0.0 
dpinf = 0.0 
machu = 0.0 
pstatu = ptot 
qu = 0.0 
rhou = 0.0 
rnftu = 0.0 
velu = 0.0 

endif 

* initial values for corrected tunnel parameters 

mach = machu 
pstat = pstatu 
q = qu 
rho = rhou 
rnft = rnftu 
vel = velu 
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Flowmeters 


Input Variables 

flodp# - flowmeter differential pressure (psi) 
flopi#- flowmeter inlet static pressure (psi) 
flotr# - flowmeter temperature (degR) 
ptot - settling chamber total pressure (psf) 
tr - ambient temperature (degR) 

Output Variables 

flbeta# - flowmeter diameter ratio 

floasq# - cross-sectional area of throat (sqft) 

floe# - discharge coefficient 

flodi# - flowmeter inlet diameter (in) 

flodt# - flowmeter throat diameter (in) 

flof# - velocity of approach factor 

flomu# - flowmeter fluid viscosity (lbsec/ft) 

flopt# - flowmeter throat static pressure (psi) 

floptn# - normalized ambient pressure (psi) 

flor# - flowmeter pressure ratio 

Horn# - flowmeter reynolds number 

flotf# - flowmeter temperature (degF) 

flottn# - normalized ambient temperature (degR) 

floy# - flowmeter expansion factor 

flsfmc# - flowmeter super compressibility factor 

wtflo# - mass flow rate (lbs/sec) 

wtflon# - normalized mass flow rate (lbs/sec) 

* note # is flowmeter number 

* maximum number of flowmeters = 15 

* the input and output variables aboved are per flowmeter 

* note The source code routine, flowmeters. F, is used to perform the flowmeters’ calculations. 
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Jet Exhaust Measurements 


Input Variables 

fml - primary flow flowmeter frequency; air system 1 (Hz) 

fm2 - primary flow flowmeter frequency; air system 2 (Hz) 

fms - secondary flow flowmeter frequency (Hz) 

machu - uncorrected free stream mach number 

pbll - tertiary flow static pressure; probe 1 (psi) 

pbl2 - tertiary flow static pressure, probe 2 (psi) 

pbl3 - tertiary flow static pressure; probe3 (psi) 

pbl4 - tertiary flow static pressure; probe4 (psi) 

pchl - plenum chamber total pressure; engine 1 (psi) 

pch2 - plenum chamber total pressure; engine 2 (psi) 

pch3 - plenum chamber total pressure; engine 3 (psi) 

pch4 - plenum chamber total pressure; engine 4 (psi) 

pfml - primary flow flowmeter pressure; air system 1 (psi) 

pfm2 - primary flow flowmeter pressure; air system 2 (psi) 

pfms - secondary flow flowmeter pressure; (psi) 

pstatu - uncorrected test section static pressure (psf) 

psv - venturi static pressure (psi) 

psi - primary jet static pressure; engine 1 (psi) 

ps2 - primary jet static pressure; engine 2 (psi) 

ps3 - primary jet static pressure; engine 3 (psi) 

ps4 - primary jet static pressure; engine 4 (psi) 

ptbll - tertiary flow total pressure; probe 1 (psi) 

ptbl2 - tertiary flow total pressure; probe 2 (psi) 

ptbl3 - tertiary flow total pressure; probe 3 (psi) 

ptbl4 - tertiary flow total pressure; probe 4 (psi) 

ptjl.l - primary jet total pressure; engine 1, probe 1 (psi) 

ptj 1.2- primary jet total pressure; engine 1 , probe 2 (psi) 

ptjl.3 - primary jet total pressure; engine 1, probe 3 (psi) 

ptj 1.4- primary jet total pressure; engine 1 , probe 4 (psi) 

ptj 1.5 - primary jet total pressure; engine 1, probe 5 (psi) 

ptj 1 .6 - primary jet total pressure; engine 1 , probe 6 (psi) 

ptj 1.7- primary jet total pressure; engine 1 , probe 7 (psi) 

ptj 1.8 - primary jet total pressure; engine 1, probe 8 (psi) 

ptj 1.9- primary jet total pressure; engine 1 , probe 9 (psi) 

ptj 1 . 10 - primary jet total pressure; engine 1 , probe 10 (psi) 

ptj 1 . 1 1 - primary jet total pressure; engine 1 , probe 1 1 (psi) 

ptj 1 . 12 - primary jet total pressure; engine 2, probe 12 (psi) 

ptj 2 . 1 - primary jet total pressure; engine 2, probe 1 (psi) 

ptj2.2 - primary jet total pressure; engine 2, probe 2 (psi) 

ptj 2 . 3 - primary jet total pressure; engine 2, probe 3 (psi) 

ptj2.4 - primary jet total pressure; engine 2, probe 4 (psi) 

ptj2.5 - primary jet total pressure; engine 2, probe 5 (psi) 

ptj2.6 - primary jet total pressure; engine 2, probe 6 (psi) 

ptj2.7 - primary jet total pressure; engine 2, probe 7 (psi) 

ptj 2 . 8 - primary jet total pressure; engine 2, probe 8 (psi) 

ptj2.9 - primary jet total pressure; engine 2, probe 9 (psi) 
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ptj 2.10 - primary jet total pressure; engine 2, probe 10 (psi) 

ptj 2.11 - primary jet total pressure; engine 2, probe 1 1 (psi) 

ptj 2.12 - primary jet total pressure; engine 2, probe 12 (psi) 

ptj 3 . 1 - primary jet total pressure; engine 3, probe 1 (psi) 
ptj3.2 - primary jet total pressure; engine 3, probe 2 (psi) 
ptj 3 . 3 - primary jet total pressure; engine 3, probe 3 (psi) 
ptj3.4 - primary jet total pressure; engine 3, probe 4 (psi) 
ptj3.5 - primary jet total pressure; engine 3, probe 5 (psi) 
ptj3.6 - primary jet total pressure; engine 3, probe 6 (psi) 
ptj3.7 - primary jet total pressure; engine 3, probe 7 (psi) 
ptj 3 . 8 - primary jet total pressure; engine 3, probe 8 (psi) 
ptj3.9 - primary jet total pressure; engine 3, probe 9 (psi) 
ptj 3. 10 - primary jet total pressure; engine 3, probe 10 (psi) 

ptj 3.11 - primary jet total pressure; engine 3, probe 1 1 (psi) 

ptj 3. 12 - primary jet total pressure; engine 3, probe 12 (psi) 

ptj4.1 - primary jet total pressure; engine 4, probe 1 (psi) 
ptj4.2 - primary jet total pressure; engine 4, probe 2 (psi) 
ptj4.3 - primary jet total pressure; engine 4, probe 3 (psi) 
ptj4.4 - primary jet total pressure; engine 4, probe 4 (psi) 
ptj4.5 - primary jet total pressure; engine 4, probe 5 (psi) 
ptj4.6 - primary jet total pressure; engine 4, probe 6 (psi) 
ptj4.7 - primary jet total pressure; engine 4, probe 7 (psi) 
ptj4.8 - primary jet total pressure; engine 4, probe 8 (psi) 
ptj4.9 - primary jet total pressure; engine 4, probe 9 (psi) 
ptj4.10 - primary jet total pressure; engine 4, probe 10 (psi) 

ptj4. 1 1 - primary jet total pressure; engine 4, probe 1 1 (psi) 

ptj4.12 - primary jet total pressure; engine 4, probe 12 (psi) 

ptot - settling chamber total pressure (psf) 
ptsl - secondary flow total pressure; probe 1 (psi) 
pts2 - secondary flow total pressure; probe 2 (psi) 
pts3 - secondary flow total pressure; probe 3 (psi) 
pts4 - secondary flow total pressure; probe 4 (psi) 
ptv - venturi total pressure (psi) 

pvenl.l - multiple critical static pressure; air system 1, upstream of 
venturi throat (psi) 

pvenl.2 - multiple critical static pressure; air system 1, downstream of 
venturi throat (psi) 

pvenl.3 - multiple critical static pressure; air systeml, upstream of 
venturi throat (psi) 

pvenl.4 - multiple critical static pressure; air systeml, downstream of 
venturi throat (psi) 

pven2.1 - multiple critical static pressure; air system 2, upstream of 
venturi throat (psi) 

pven2.2 - multiple critical static pressure; air system 2, downstream of 
venturi throat (psi) 

pven2.3 - multiple critical static pressure; air system 2, upstream of 
venturi throat (psi) 

pven2.4 - multiple critical static pressure; air system 2, downstream of 
venturi throat (psi) 

pvril.l - in-line(not mcv) venturi static pressure; air system 1, venturi 1 (psi) 
pvril.2 - in-line(not mcv) venturi static pressure; air system 1, venturi 2 (psi) 

13 


pvril.3 - in-line(not mcv) venturi static pressure; air system 1, venturi 3 (psi) 

pvril.4 - in-line(not mcv) venturi static pressure; air system 1, venturi 4 (psi) 

pvri2.1 - in-line(not mcv) venturi static pressure; air system 2, venturi 1 (psi) 

pvri2.2 - in-line(not mcv) venturi static pressure; air system 2, venturi 2 (psi) 

pvri2.3 - in-line(not mcv) venturi static pressure; air system 2, venturi 3 (psi) 

pvri2.4 - in-line(not mcv) venturi static pressure; air system 2, venturi 4 (psi) 

qu - uncorrected tunnel dynamic pressure with compressibility (psf) 

tchl - plenum chamber total temperature; engine 1 (degF) 

tch2 - plenum chamber total temperature; engine 2 (degF) 

tch3 - plenum chamber total temperature; engine 3 (degF) 

tch4 - plenum chamber total temperature; engine 4 (degF) 

tfml - primary flow flowmeter temperature; air system 1 (degF) 

tfm2 - primary flow flowmeter temperature; air system 2 (degF) 

tfms - secondary flow flowmeter temperature (degF) 

ttbl - tertiary flow total temperature (degF) 

ttjl.l - primary jet total temperature; engine 1; probe 1 (degF) 

ttjl.2 - primary jet total temperature; engine 1; probe 2 (degF) 

ttjl.3 - primary jet total temperature; engine 1; probe 3 (degF) 

ttj 1.4- primary jet total temperature; engine 1 ; probe 4 (degF) 

ttjl.5 - primary jet total temperature; engine 1, probe 5 (degF) 

ttj 1.6 - primary jet total temperature; engine 1, probe 6 (degF) 

ttj2. 1 - primary jet total temperature, engine 2, probe 1 (degF) 

ttj2.2 - primary jet total temperature; engine 2, probe 2 (degF) 

ttj2.3 - primary jet total temperature; engine 2, probe 3 (degF) 

ttj2.4 - primary jet total temperature; engine 2, probe 4 (degF) 

ttj2.5 - primary jet total temperature; engine 2, probe 5 (degF) 

ttj2.6 - primary jet total temperature; engine 2, probe 6 (degF) 

ttj3.1 - primary jet total temperature; engine 3, probe 1 (degF) 

ttj3.2 - primary jet total temperature; engine 3, probe 2 (degF) 

ttj3.3 - primary jet total temperature; engine 3, probe 3 (degF) 

ttj3.4 - primary jet total temperature; engine 3, probe 4 (degF) 

ttj 3 . 5 - primary jet total temperature; engine 3, probe 5 (degF) 

ttj3.6 - primary jet total temperature; engine 3, probe 6 (degF) 

ttj4. 1 - primary jet total temperature; engine 4, probe 1 (degF) 

ttj4.2 - primary jet total temperature; engine 4, probe 2 (degF) 

ttj4.3 - primary jet total temperature; engine 4, probe 3 (degF) 

ttj4.4 - primary jet total temperature; engine 4, probe 4 (degF) 

ttj4.5 - primary jet total temperature; engine 4, probe 5 (degF) 

ttj4.6 - primary jet total temperature; engine 4, probe 6 (degF) 

ttsec - secondary flow total temperature; (degF) 

ttv - tertiary flow venturi temperature (degF) 

tv 1.1 - venturi temperature; air system 1, venturi 1 (degF) 

tv 1.2 - venturi temperature; air system 1, venturi 2 (degF) 

tv2.1 - venturi temperature; air system 2, venturi 1 (degF) 

tv2.2 - venturi temperature; air system 2, venturi 2 (degF) 

tvril.l - in-line(not mcv) venturi temperature; air system 1, venturi 1 (degF) 

tvril.2 - in-line(not mcv) venturi temperature; air system 1, venturi 2 (degF) 

tvril.3 - in-line(not mcv) venturi temperature; air system 1, venturi 3 (degF) 

tvril.4 - in-line(not mcv) venturi temperature; air system 1, venturi 4 (degF) 

tvri2.1 - in-line(not mcv) venturi temperature; air system 2, venturi 1 (degF) 

tvri2.2 - in-line(not mcv) venturi temperature; air system 2, venturi 2 (degF) 
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tvri2.3 - in-line(not mcv) venturi temperature; air system 2, venturi 3 (degF) 
tvri2.4 - in-line(not mcv) venturi temperature; air system 2, venturi 4 (degF) 

Output Variables 

anlzl - total throat area; air system 1 
anlz2 - total throat area; air system 2 

cfi 1 - ideal thrust coefficient based on mass flow rate measured by flowmeter; 
air system 1 

cfi2 - ideal thrust coefficient based on mass flow rate measured by flowmeter; 
air system 2 

cfichrl - ideal thrust coefficient based on mass flow rate computed from plenum 
chamber measurements; air system 1 

cfichr2 - ideal thrust coefficient based on mass flow rate computed from plenum 
chamber measurements; air system 2 
fi 1 - ideal thrust of total primary exhaust system based on mass flow rate 
measured by flowmeter; air system 1 (lbs) 
fi2 - ideal thrust of total primary exhaust system based on mass flow rate 
measured by flowmeter; air system 2 (lbs) 
fia - average ideal thrust of total primary exhaust system (lbs) 
fichrl - ideal thrust of total primary exhaust system based on mass flow rate 
computed from plenum chamber measurements; air system 1 (lbs) 
fichr2 - ideal thrust of total primary exhaust system based on mass flow rate 
computed from plenum chamber measurements; air system 2 (lbs) 
fiengl - ideal thrust based on mass flow rate computed from individual 
plenum chamber measurements; engine 1 (lbs) 
fieng2 - ideal thrust based on mass flow rate computed from individual 
plenum chamber measurements; engine 2 (lbs) 
fieng3 - ideal thrust based on mass flow rate computed from individual 
plenum chamber measurements; engine 3 (lbs) 
fieng4 - ideal thrust based on mass flow rate computed from individual 
plenum chamber measurements; engine 4 (lbs) 
mbldot - tertiary flow mass flow rate (slugs/sec) 
mdot 1 - primary flow mass flow rate measured by flowmeter; 
air system 1 (slugs/sec) 

mdot2 - primary flow mass flow rate measured by flowmeter; 
air system 2 (slugs/sec) 

mdota - average primary flow mass flow rate (slugs/sec) 
mdotchl - primary flow mass flow rate computed from plenum chamber 
measurements; air system 1 (slugs/sec) 
mdotch 2 - primary flow mass flow rate computed from plenum chamber 
measurements; air system 2 (slugs/sec) 
mductl - mach number; engine 1 
mduct2 - mach number; engine 2 
mduct3 - mach number; engine 3 
mduct4 - mach number; engine 4 
msdot - secondary flow mass flow rate (slugs/sec) 
pblave - average tertiary flow static pressure (psi) 
pchoke - primary choked flow jet total pressure ratio 
pjetl - primary jet pressure; air system 1 (psi) 
pjet2 - primary jet pressure; air system 2 (psi) 
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psec - average secondary flow static pressure (psi) 
ptblav - average tertiary flow total pressure (psi) 
ptb/pj - ratio of tertiary flow total pressure to primary jet total pressure 
ptb/po - ratio of tertiary flow total pressure to free stream total pressure 
ptengl - average primary jet total pressure; engine 1 (psi) 
pteng2 - average primary jet total pressure; engine 2 (psi) 
pteng3 - average primary jet total pressure; engine 3 (psi) 
pteng4 - average primary jet total pressure; engine 4 (psi) 
ptengol - ratio of average primary jet total pressure in engine 1 
to tunnel static pressure 

ptengo2 - ratio of average primary jet total pressure in engine 2 
to tunnel static pressure 

ptengo3 - ratio of average primary jet total pressure in engine 3 
to tunnel static pressure 

ptengo4 - ratio of average primary jet total pressure in engine 4 
to tunnel static pressure 

ptj/pol - primary jet total pressure ratio (all engines); air system 1 
ptj/po2 - primary jet total pressure ratio (all engines); air system 2 
ptj/poa - average primary jet total pressure ratio (all engines) 
ptsec - average secondary flow total pressure (psi) 

pts/po - ratio of secondary flow total pressure to free stream total pressure 

pts/ptj - ratio of secondary flow total pressure to primary jet total pressure 

pval - coefficient pvl; air system 1 

pva2 - coefficient pvl; air system 2 

tasl - temperature; air system 1 

tas2 - temperature; air system 2 

thetase - secondary flow corrected mass flow ratio (slugs/sec) 
thetbl - tertiary flow corrected mass flow rate (slugs/sec) 
ttengl - average primary jet total temperature; engine 1 (degF) 
tteng2 - average primary jet total temperature; engine 2 (degF) 
tteng3 - average primary jet total temperature; engine 3 (degF) 
tteng4 - average primary jet total temperature; engine 4 (degF) 
ttjav - average primary jet total temperature; all engines and air systems (degF) 
ttjavgl - average primary jet total temperature; all engines, 
air system 1 (degF) 

ttjavg2 - average primary jet total temperature; all engines, 
air system 2 (degF) 

tval - average venturi temperature; air system 1 (degF) 
tva2 - average venturi temperature; air system 2 (degF) 
vratiol - ratio of multiple critical venturi static pressures; air system 1; 
should be less than 0.93 

vratio2 - ratio of multiple critical venturi static pressures; air system 2; 
should be less than 0.93 

wil - ideal weight flow rate; air system 1 (lbs/sec) 

wi2 - ideal weight flow rate; air system 2 (lbs/sec) 

wiengl - ideal weight flow rate; engine 1 (lbs/sec) 

wieng2 - ideal weight flow rate; engine 2 (lbs/sec) 

wieng3 - ideal weight flow rate; engine 3 (lbs/sec) 

wieng4 - ideal weight flow rate; engine 4 (lbs/sec) 

wmcvl - multiple critical venturi weight flow rate; air system 1 (lbs/sec) 

wmcv2 - multiple critical venturi weight flow rate; air system 2 (lbs/sec) 
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wmcv/wil - ratio of multiple critical venturi weight flow rate to ideal weight 
flow rate; air system 1 

wmcv/wi2 - ratio of multiple critical venturi weight flow rate to ideal weight 
flow rate; air system 2 

wpl - primary flow flowmeter weight flow rate; air system 1 (lbs/sec) 
wp2 - primary flow flowmeter weight flow rate; air system 2 (lbs/sec) 
wpbl - tertiary flow venturi weight flow rate (lbs/sec) 

wpchrl - total primary flow weight flow rate calculated from plenum chamber 
measurements; air system 1 (lbs/sec) 

wpchr2 - total primary flow weight flow rate calculated from plenum chamber 
measurements; air system 2 (lbs/sec) 
wpch/wil - total primary flow discharge coefficient computed from 
plenum chamber measurements; air system 1 
wpch/wi2 - total primary flow discharge coefficient computed from 
plenum chamber measurements; air system 2 
wpengl - primary flow weight flow rate computed from plenum chamber 
measurements; engine 1 (lbs/sec) 

wpeng2 - primary flow weight flow rate computed from plenum chamber 
measurements; engine 2 (lbs/sec) 

wpeng3 - primary flow weight flow rate computed from plenum chamber 
measurements; engine 3 (lbs/sec) 

wpeng4 - primary flow weight flow rate computed from plenum chamber 
measurements; engine 4 (lbs/sec) 
wpsec - secondary flow weight flow rate (lbs/sec) 

wp/wil - primary flow discharge coefficient using flowmeter weight flow rate; 
air system 1 

wp/wi2 - primary flow discharge coefficient using flowmeter weight flow rate; 
air system 2 

wp/wiel - discharge coefficient computed from plenum chamber 
measurements; engine 1 

wp/wie2 - discharge coefficient computed from plenum chamber 
measurements; engine 2 

wp/wie3 - discharge coefficient computed from plenum chamber 
measurements; engine 3 

wp/wie4 - discharge coefficient computed from plenum chamber 
measurements; engine 4 

wvril.l - in-line venturi weight flow rate; air system 1, venturi 1 (lbs/sec) 
wvril.2 - in-line venturi weight flow rate; air system 1, venturi 2 (lbs/sec) 
wvril.3 - in-line venturi weight flow rate; air system 1, venturi 3 (lbs/sec) 
wvril.4 - in-line venturi weight flow rate; air system 1, venturi 4 (lbs/sec) 
wvri2. 1 - in-line venturi weight flow rate; air system 2, venturi 1 (lbs/sec) 
wvri2.2 - in-line venturi weight flow rate; air system 2, venturi 2 (lbs/sec) 
wvri3.3 - in-line venturi weight flow rate; air system 2, venturi 3 (lbs/sec) 
wvri4.4 - in-line venturi weight flow rate; air system 2, venturi 4 (lbs/sec) 


* note The source code routine, moduli. F, is used to perform the jet exhaust 

* measurements ’ calculations. 
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Refer to pages 14 thru 43 of reference 5 for a detailed explanation of the jet exhaust calculations. 
Also, refer to reference 1 for information concerning the Multiple Critical Venturi System. 
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Balance Loads and Model Attitudes 


First, calculations are performed for initial loads and second order interactions due to initial loads. 
Next, corrections for first and second order balance interactions, high interactions/high model restraints, 
air line pressure, weight tares, method of attachment, and sting deflections are performed. Axis rotations 
follow these corrections. 

Axis rotations include gravity axis to balance axis, gravity axis to model(body) axis, wind axis to 
gravity axis, and wind axis to model(body) axis. Angle calculations are performed for the balance and 
model(body) axes. 

Next, calculations for the model height, balance components rotated and translated to the 
model(body) axis, applying estimates of Heyson’s boundary interference, blockage and jet boundary 
corrections, Fleyson’s boundary(wall) interference, and base(cavity) pressures are performed. Finally, 
balance components and coefficients are calculated for the model, stability, wind, and reference axes. 

Although only calculations for the balance one loads are explained, a maximum of four balances 
may be used. Unless noted, the calculations for balances two thru four are the same as the calculations for 
balance one. The variable names used in the calculations for each balance are named according to the 
balance number. 
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Initial Loads 


Input Variables 

phiO - model roll angle at wind-off zero (deg) 

thetaO - model pitch angle at wind-off zero (deg) 

wafl - axial force attitude load (lbs) 

wnf 1 - axial force attitude load (lbs) 

wsfl - side force attitude load (lbs) 

wxpml - pitching moment attitude load (in-lbs) 

wxyml - yawing moment attitude load (in-lbs) 

wyrml - rolling moment attitude load (in-lbs) 

wyyml - yawing moment attitude load (in-lbs) 

wzpml - pitching moment attitude load (in-lbs) 

wzrml - rolling moment attitude load (in-lbs) 

Output Variables 

afl.O - axial force initial load (lbs) 

nf 1 .0 - normal force initial load (lbs) 

pml.O - pitching moment initial load (in-lbs) 

rml.O - rolling moment initial load (in-lbs) 

sf 1 .0 - side force initial load (lbs) 

yml.O - yawing moment initial load (in-lbs) 


* note dtor = 0.0174532925199 (converts degrees to radians} 

afl.O = wafl * sin(thetaO*dtor) 
sfl.O = wsfl * sin(phiO*dtor) * cos(thetaO*dtor) 
nfl.O = -wnfl * cos(phiO*dtor) * cos(thetaO*dtor) 
rml.O = (wyrml * cos(phiO*dtor) * cos(thetaO*dtor)) + 

(wzrml * sin(phiO*dtor) * cos(thetaO*dtor)) 
pml.O = (-wxpml * cos(phiO*dtor) * cos(thetaO*dtor)) + 
(wzpml * sin(thetaO*dtor)) 

yml.O = (wxyml * sin(phiO*dtor) * cos(thetaO*dtor)) + 
(wyyml * sin(thetaO*dtor)) 
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Second Order Interactions Due To Initial Loads 


Input Variables 

afl.O - axial force initial load (lbs) 

nfl.O - normal force initial load (lbs) 

pm 1.0 - pitching moment initial load (in-lbs) 

rml.O - rolling moment initial load (in-lbs) 

sf 1 .0 - side force initial load (lbs) 

yml.O - yawing moment initial load (in-lbs) 

Output Variables 

afezl - axial force second order interactions due to initial loads (lbs) 
nfezl- normal force second order interactions due to initial loads (lbs) 
pmezl- pitching moment second order interactions due to initial loads (in-lbs) 
rmezl- rolling moment second order interactions due to initial loads (in-lbs) 
sfez 1 - side force second order interactions due to initial loads (lbs) 
ymezl- yawing moment second order interactions due to initial loads (in-lbs) 


* note execute subfunction ctrnl which computes the second order interactions 

* due to initial loads. 

* note afezl, sfezl, nfezl, rmezl, pmezl, ymezl are not written to output datafile 
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First & Second Order Balance Interactions 


Input Variables 

afezl - axial force second order interactions due to initial loads(lbs) 
af 1 - uncorrected axial force (lbs) 
afl.O - axial force initial load (lbs) 

nfezl- normal force second order interactions due to initial loads(lbs) 
nf 1 - uncorrected normal force (lbs) 
nfl.O - normal force initial load (lbs) 

pmezl- pitching moment second order interactions due to initial loads(in-lbs) 
pml - uncorrected pitching moment (in-lbs) 
pm 1.0 - pitching moment initial load (in-lbs) 

rmezl- rolling moment second order interactions due to initial loads(in-lbs) 
rml - uncorrected rolling moment (in-lbs) 
rml.O - rolling moment initial load (in-lbs) 

sfezl- side force second order interactions due to initial loads(lbs) 
sf 1 - uncorrected side force (lbs) 
sf 1 .0 - side force initial load (lbs) 

ymezl- yawing moment second order interactions due to initial loads(in-lbs) 
yml - uncorrected yawing moment (in-lbs) 
yml.O - yawing moment initial load (in-lbs) 

Output Variables 

afl.2 - axial force load corrected for 1 st & 2 nd order balance interactions (lbs) 

nfl.2 - normal force load corrected for 1 st & 2 nd order balance interactions (lbs) 

pm 1.2 - pitching moment load corrected for 1 st & 2 nd order balance interactions (in-lbs) 

rml. 2 - rolling moment load corrected for 1 st & 2 nd order balance interactions (in-lbs) 

sfl.2 - side force load corrected for 1 st & 2 nd order balance interactions (lbs) 

yml. 2 - yawing moment load corrected for 1 st & 2 nd order balance interactions (in-lbs) 


* note execute subfunction cintr which computes corrected delta loads for first 

* and second order balance interactions. 
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High Interactions and High Model Restraints 


The six balance loads are corrected for high interactions and high model restraints. A 6x6 balance 

interaction matrix is used for the interaction factors. 

Input Variables 

afl.2 - axial force load corrected for 1 st & 2 nd order balance interactions (lbs) 

blaal - axial force(axial) interaction factor; constant; default value is 1.0 

blnal - axial force(normal) interaction factor; constant; default value is 0.0 

blpal - axial force(pitch) interaction factor; constant; default value is 0.0 

blral - axial force(roll) interaction factor; constant; default value is 0.0 

blsal - axial force(side) interaction factor; constant; default value is 0.0 

blyal - axial force(yaw) interaction factor; constant; default value is 0.0 

blanl - normal force(axial) interaction factor; constant; default value is 0.0 

blnnl - normal force(normal) interaction factor; constant; default value is 1.0 

blpnl - normal force(pitch) interaction factor; constant; default value is 0.0 

blrnl - normal force(roll) interaction factor; constant; default value is 0.0 

blsnl - normal force(side) interaction factor; constant; default value is 0.0 

blynl - normal force(yaw) interaction factor; constant; default value is 0.0 

blapl - pitching moment(axial) interaction factor; constant; default value is 0.0 

blnpl - pitching moment(normal) interaction factor; constant; default value is 0.0 

blppl - pitching moment(pitch) interaction factor; constant; default value is 1.0 

blrpl - pitching moment(roll) interaction factor; constant; default value is 0.0 

blspl - pitching moment(side) interaction factor; constant; default value is 0.0 

blypl - pitching moment(yaw) interaction factor; constant; default value is 0.0 

b lari - rolling moment(axial) interaction factor; constant; default value is 0.0 

blmT - rolling moment(normal) interaction factor; constant; default value is 0.0 

blprl - rolling moment(pitch) interaction factor; constant; default value is 0.0 

blrrl - rolling moment(roll) interaction factor; constant; default value is 1.0 

blsrl - rolling moment(side) interaction factor; constant; default value is 0.0 

blyrl - rolling moment(yaw) interaction factor; constant; default value is 0.0 

b lasl - side force(axial) interaction factor; constant; default value is 0.0 

blnsl - side force(normal) interaction factor; constant; default value is 0.0 

blpsl - side force(pitch) interaction factor; constant; default value is 0.0 

blrsl - side force(roll) interaction factor; constant; default value is 0.0 

blssl - side force(side) interaction factor; constant; default value is 1.0 

blysl - side force(yaw) interaction factor; constant; default value is 0.0 

blayl - yawing moment(axial) interaction factor; constant; default value is 0.0 

blnyl - yawing moment(normal) interaction factor; constant; default value is 0.0 

blpyl - yawing moment(pitch) interaction factor; constant; default value is 0.0 

blryl - yawing moment(roll) interaction factor; constant; default value is 0.0 

blsyl - yawing moment(side) interaction factor; constant; default value is 0.0 

blyyl - yawing moment(yaw) interaction factor; constant; default value is 1.0 

nfl.2 - normal force load corrected for 1 st & 2 nd order balance interactions (lbs) 

pm 1.2 - pitching moment load corrected for 1 st & 2 nd order balance interactions (in-lbs) 

rml.2 - rolling moment load corrected for 1 st & 2 nd order balance interactions (in-lbs) 

sfl.2 - side force load corrected for 1 st & 2 nd order balance interactions (lbs) 

yml.2 - yawing moment load corrected for 1 st & 2 nd order balance interactions (in-lbs) 
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Output Variables 


afl.3 - axial force load corrected for high interactions and high model retraints (lbs) 
nfl.3 - normal force load corrected for high interactions and high model restraints (lbs) 
pm 1.3 - pitching moment load corrected for high interactions and high model restraints (in-lbs) 
rml.3 - rolling moment load corrected for high interactions and high model restraints (in-lbs) 
sfl.3 - side force load corrected for high interactions and high model restraints (lbs) 
yml.3 - yawing moment load corrected for high interactions and high model restraints (in-lbs) 


afl.3 = (afl.2*blaal) + (sfl.2*blsal) + (nfl.2*blnal) +(rml.2*blral) + (pml.2*blpal) 

+ (yml.2*blyal) 

sfl.3 = (afl.2*blasl) + (sfl.2*blssl) + (nfl.2*blnsl) + (rml.2*blrsl) + (pml.2*blpsl) 

+ (yml.2*blysl) 

nfl.3 = (afl.2*blanl) + (sfl.2*blsnl) + (nfl.2*blnnl) + (rml.2*blrnl) + (pml.2*blpnl) 
+ (yml.2*blynl) 

rml.3 = (afl.2*blarl) + (sfl.2*blsrl) + (nfl.2*blnrl) + (rml.2*blrrl) + (pml.2*blprl) 

+ (yml.2*blyrl) 

pml.3 = (afl.2*blapl) + (sfl.2*blspl) + (nfl.2*blnpl) + (rml.2*blrpl) + (pml.2*blppl) 
+ (yml.2*blypl) 

yml.3 = (afl.2*blayl) + (sfl.2*blsyl) + (nfl.2*blnyl) + (rml.2*blryl) + (pml.2*blpyl) 
+ (yml.2*blyyl) 
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Air Line Pressure Corrections 


The six balance loads are corrected for air line pressure. 

Input Variables 

afl.3 - axial force load corrected for high interactions and high model retraints (lbs) 

nfl.3 - normal force load corrected for high interactions and high model restraints (lbs) 

pcal.l - axial force air line pressure coefficient 1; constant; default value is 0.0 

pca2.1 - axial force air line pressure coefficient 2; constant; default value is 0.0 

pca3.1 - axial force air line pressure coefficient 3; constant; default value is 0.0 

pen 1.1 - normal force air line pressure coefficient 1; constant; default value is 0.0 

pcn2.1 - normal force air line pressure coefficient 2; constant; default value is 0.0 

pcn3. 1 - normal force air line pressure coefficient 3; constant; default value is 0.0 

pep 1.1 - pitching moment air line pressure coefficient 1; constant; default value is 0.0 

pcp2.1 - pitching moment air line pressure coefficient 2; constant; default value is 0.0 

pcp3. 1 - pitching moment air line pressure coefficient 3; constant; default value is 0.0 

pcrl.l - rolling moment air line pressure coefficient 1; constant; default value is 0.0 

pcr2.1 - rolling moment air line pressure coefficient 2; constant; default value is 0.0 

pcr3.1 - rolling moment air line pressure coefficient 3; constant; default value is 0.0 

pcs 1.1 - side force air line pressure coefficient 1; constant; default value is 0.0 

pcs2. 1 - side force air line pressure coefficient 2; constant; default value is 0.0 

pcs3.1 - side force air line pressure coefficient 3; constant; default value is 0.0 

pcyl.l - yawing moment air line pressure coefficient 1; constant; default value is 0.0 

pcy2.1 - yawing moment air line pressure coefficient 2; constant; default value is 0.0 

pcy3.1 - yawing moment air line pressure coefficient 3; constant; default value is 0.0 

pm 1.3 - pitching moment load corrected for high interactions and high model restraints (in-lbs) 

pst - air line pressure (psi) 

rml.3 - rolling moment load corrected for high interactions and high model restraints (in-lbs) 
sfl.3 - side force load corrected for high interactions and high model restraints (lbs) 
yml.3 - yawing moment load corrected for high interactions and high model restraints (in-lbs) 

Output Variables 

afl.3 - axial force load corrected for high interactions, high model retraints and air line pressure (lbs) 
nfl.3 - normal force load corrected for high interactions, high model restraints and air line 
pressure (lbs) 

pm 1.3 - pitching moment load corrected for high interactions, high model restraints, and air line 
pressure (in-lbs) 

rml.3 - rolling moment load corrected for high interactions, high model restraints, and air line 
pressure (in-lbs) 

sfl.3 - side force load corrected for high interactions, high model restraints and air line pressure (lbs) 
yml.3 - yawing moment load corrected for high interactions, high model restraints, and air line 
presssure (in-lbs) 
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afl.3 = afl.3 - ((pst*pcal.l) + (pst*pst*pca2.1) + (pst*pst*pst*pca3.1)) 
sfl.3 = sfl.3 - ((pst*pcsl.l) + (pst*pst*pcs2.1) + (pst*pst*pst*pcs3.1)) 
nfl.3 = nfl.3 - ((pst*pcnl.l) + (pst*pst*pcn2.1) +(pst*pst*pst*pcn3.1)) 
rml.3 = rml.3 - ((pst*pcrl.l) + (pst*pst*pcr2.1) +(pst*pst*pst*pcr3.1)) 
pml.3 = pml.3 - ((pst*pcpl.l) + (pst*pst*pcp2.1) +(pst*pst*pst*pcp3.1)) 
yml.3 = yml.3 - ((pst*pcyl.l) + (pst*pst*pcy2.1) +(pst*pst*pst*pcy3.1)) 
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Sting Deflections 

Input Variables 

afl.O - axial force initial load (lbs) 

afl.3 - axial force load corrected for high interactions, high model retraints 
and air line pressure (lbs) 
kdfl - sting deflection corrections flag; constant 
nfl.O - normal force initial load (lbs) 

nfl.3 - normal force load corrected for high interactions, high model restraints 
and air line pressure (lbs) 
phidal - axial force roll deflection; constant 
phidnl - normal force roll deflection; constant 
phidpl - pitching moment roll deflection; constant 
phidrl - rolling moment roll deflection; constant 
phidsl - side force roll deflection; constant 
phidyl - yawing moment roll deflection; constant 
pm 1.0 - pitching moment initial load (in-lbs) 

pm 1.3 - pitching moment load corrected for high interactions, high model 
restraints, and air line pressure (in-lbs) 
psidal -axial force yaw deflection; constant 
psidnl - normal force yaw deflection; constant 
psidpl - pitching moment yaw deflection; constant 
psidrl - rolling moment yaw deflection; constant 
psidsl - side force yaw deflection; constant 
psidyl - yawing moment yaw deflection; constant 
rml.O - rolling moment intial load (in-lbs) 

rml.3 - rolling moment load corrected for high interactions, high model 
restraints, and air line pressure (in-lbs) 
sfl.O - side force initial load (lbs) 

sfl.3 - side force load corrected for high interactions, high model restraints 
and air line pressure (lbs) 
thedal - axial force pitch deflection; constant 
thednl - normal force pitch deflection; constant 
thedpl - pitching moment pitch deflection; constant 
thedrl - rolling moment pitch deflection; constant 
thedsl - side force pitch deflection; constant 
thedyl - yawing moment pitch deflection; constant 
yml.O - yawing moment initial load (in-lbs) 

yml.3 - yawing moment load corrected for high interactions, high model 
restraints, and air line presssure (in-lbs) 

Output Variables 

phid - roll deflection angle (deg) 
psid - yaw deflection angle (deg) 
thetad - pitch deflection angle (deg) 
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note 


The default value is 0.0 for phidal, phidsl, phidnl, phidrl, phidpl, 
phidyl, thedal, thedsl, thednl, thedrl, thedpl, thedyl, psidal, psidsl, 
psidnl, psidrl, psidpl, psidyl 


* 

* 

* 


if ( kdfl = -1 ) then 

*no sting deflection corrections 

phid = 0.0 

thetad = 0.0 

psid = 0.0 

else if ( kdfl = 0 ) then 

*sting deflection corrections using total loads 

phid = (phidal *(afl .0+af 1 .3)) + (phidsl*(sfl.0+sfl.3)) + 
(phidnl*(nfl.0+nfl.3)) + (phidrl*(rml.0+rml.3)) + 
(phidpl*(pml.0+pml.3)) + (phidyl*(yml.0+yml.3)) 

thetad = (thedal*(afl.0+afl.3)) + (thedsl*(sfl.0+sfl.3)) + 
(thednl*(nfl.0+nfl.3)) + (thedrl*(rml.0+rml.3)) + 
(thedpl*(pml.0+pml.3)) + (thedyl*(yml.0+yml.3)) 

psid = (psidal*(afl.0+afl.3)) + (psids 1 *(sf 1 .0+sf 1 .3)) + 
(psidnl *(nfl.0+nf 1.3)) + (psidrl*(rml.0+rml.3)) + 
(psidpl*(pml.0+pml.3)) + (psidyl*(yml.0+yml.3)) 

else if ( kdfl = 1 ) then 

*sting deflection corrections using delta loads 

phid = (phidal*afl.3) + (phids 1 *sf 1 .3) + (phidnl*nfl.3) + 
(phidrl*rml.3) + (phidpl*pml.3) + (phidyl*yml.3) 

thetad = (thedal*afl.3) + (thedsl*sfl.3) + thednl*nfl.3) + 

(thedrl*rml.3) + (thedpl*pml.3) + (thedyl*yml.3) 

psid = (psidal*afl.3) + (psidsl*sfl.3) + (psidnl*nfl.3) + 
(psidrl*rml.3) + (psidpl*pml.3) + (psidyl*yml.3) 


endif 
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Weight Tares 

Input Variables 

afl.3 - axial force load corrected for high interactions, high model restraints 
and air line pressure (lbs) 

nfl.3 - normal force load corrected for high interactions, high model restraints 
and air line pressure (lbs) 
phi - model roll angle (deg) 
phiO - model roll angle at wind-off zero (deg) 

pm 1.3 - pitching moment load corrected for high interactions, high model 
restraints, and air line pressure (in-lbs) 
rml.3 - rolling moment load corrected for high interactions, high model 
restraints, and air line pressure (in-lbs) 
sfl.3 - side force load corrected for high interactions, high model restraints 
and air line pressure (lbs) 
theta - model pitch angle (deg) 
thetaO - model pitch angle at wind-off zero (deg) 
wafl - axial force attitude load (lbs) 
wnf 1 - axial force attitude load (lbs) 
wsfl - side force attitude load (lbs) 
wxpml - pitching moment attitude load (in-lbs) 
wxyml - yawing moment attitude load (in-lbs) 
wyrml - rolling moment attitude load (in-lbs) 
wyyml - yawing moment attitude load (in-lbs) 
wzpml - pitching moment attitude load (in-lbs) 
wzrml - rolling moment attitude load (in-lbs) 

yml.3 - yawing moment load corrected for high interactions, high model 
restraints, and air line presssure (in-lbs) 

Output Variables 

aftare - axial force tare correction (lbs) 

afl.4 - axial force load corrected for weight tares (lbs) 

nftare - normal force tare correction (lbs) 

nfl.4 - normal force load corrected for weight tares (lbs) 

pmtare - pitching moment tare correction (in-lbs) 

pm 1.4 - pitching moment load corrected for weight tares (in-lbs) 

rmtare - rolling moment tare correction (in-lbs) 

rml.4 - rolling moment load corrected for weight tares (in-lbs) 

sftare - side force tare correction (lbs) 

sfl.4 - side force load corrected for weight tares (lbs) 

ymtare - yawing moment tare correction (in-lbs) 

yml.4 - yawing moment load corrected for weight tares (in-lbs) 
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* note 


dtor = 0.0174532925199 (converts degrees to radians} 


aftare = wafl * (sin(theta*dtor) - sin(thetaO*dtor)) 
sftare = wsfl * (cos(theta*dtor) * sin(phi*dtor)) - 
(cos(thetaO*dtor) * sin(phiO*dtor)) 
nftare = -wnfl* (cos(theta*dtor) * cos(phi*dtor)) - 
(cos(thetaO*dtor) * cos(phiO*dtor)) 

if (wafl = 0.0 & wsfl = 0.0 & wnfl = 0.0) then 
rmtare = 0.0 
pmtare - 0.0 
ymtare = 0.0 
else 

rmtare = ((wzrml/wsfl)*sftare) - ((wyrml/wnfl)*nftare) 
pmtare = ((wzpml/wafl)*aftare) + ((wxpml/wnfl)*nftare) 
ymtare = ((wxyml/wsfl)*sftare) + ((wyyml/wafl)*aftare) 
endif 

afl.4 = afl.3 - aftare 
sfl.4 = sfl.3 - sftare 
nf 1 .4 = nf 1 .3 - nftare 
rml.4 = rmfl.3 - rmtare 
pml.4 = pmfl.3 - pmtare 
yml.4 = yml.3 - ymtare 
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Method of Attachment 


Input Variables 

afl.4 - axial force load corrected for weight tares (lbs) 

ksignl - balance attachment; constant ( 1 = normal balance attachment, 

-1 = grounding balance by opposite end) 
nfl.4 - normal force load corrected for weight tares (lbs) 
pm 1.4 - pitching moment load corrected for weight tares (in-lbs) 
rml.4 - rolling moment load corrected for weight tares (in-lbs) 
sfl.4 - side force load corrected for weight tares (lbs) 
yml.4 - yawing moment load corrected for weight tares (in-lbs) 

Output Variables 

afl.5 - axial force load corrected for method of attachment (lbs) 

nfl.5 - normal force load corrected for method of attachment (lbs) 

pm 1.5 - pitching moment load corrected for method of attachment (in-lbs) 

rml.5 - rolling moment load corrected for method of attachment (in-lbs) 

sfl.5 - side force load corrected for method of attachment (lbs) 

yml.5 - yawing moment load corrected for method of attachment (in-lbs) 


afl.5 = afl.4 * ksignl 
sfl.5 = sfl.4 * ksignl 
nfl.5 = nfl.4 * ksignl 
rml.5 = rml.4 * ksignl 
pml.5 = pml.4 * ksignl 
yml.5 = yml.4 * ksignl 
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Angles’ Calculations 


Subfunction euler makes use of the angles’ rotation flags and the 3 by 3 identity, yaw, pitch, and 
roll matrices to rotate(transform) angles from one axis to another axis. Euler is used to calculate the 
rotation matrices for the gravity axis to balance axis, gravity axis to model(body) axis, and wind axis 
to gravity axis rotations. These matrices are used to generate the angles’ calculations. 


Identity matrix(3x3) 


1.0 

0.0 

0.0 

0.0 

1.0 

0.0 

0.0 

0.0 

1.0 


Yaw matrix(3x3) 


cos(angle) 

-sin(angle) 

0.0 

sin( angle) 

cos(angle) 

0.0 

0.0 

0.0 

1.0 


Pitch matrix(3x3) 


cos(angle) 

0.0 

-sin(angle) 

0.0 

1.0 

0.0 

sin( angle) 

0.0 

cos(angle) 

Roll matrix(3x3) 


1.0 

0.0 

0.0 

0.0 

cos(angle) 

-sin(angle) 

0.0 

sin(angle) 

cos(angle) 
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Rotation of Angles from Gravity Axis to Balance Axis 

Input Variables 

gbangles - nine-element array of angles (Euler input) 

gbflg - nine-element array of rotation flags (Euler input); one flag for each angle 
where values are [ - 1 = yaw rotation, 0 = pitch rotation, 1 = roll rotation ] 
idmat - 3x3 identity array (Euler input) 
phid - deflection roll angle (deg) 
phik - knuckle roll angle (deg) 
phis - strut roll angle (deg) 
psid - deflection yaw angle (deg) 
psik - knuckle yaw angle (deg) 
psis - strut yaw angle (deg) 
thetad - deflection pitch angle (deg) 
thetak - knuckle pitch angle (deg) 
thetas - strut pitch angle (deg) 

Output Variables 

gbmat - 3x3 gravity to balance axis array (Euler output) 


* note idmat and gbmat are double precision 

* note idmat, gbangles, gbflg, gbmat are not written to output datafile 

For a yaw-pitch-roll rotation, Let 

gbangles(l) = psis 
gbangles(2) = thetas 
gbangles(3) = phis 
gbangles(4) = psik 
gbangles(5) = thetak 
gbangles(6) = phik 
gbangles(7) = psid 
gbangles(8) = thetad 
gbangles(9) = phid 
gbflg(l) = -1 
gbflg(2) = 0 
gbflg(3) = 1 
gbflg(4) = -1 
gbflg(5) = 0 
gbflg(6) = 1 
gbflg(7) = -1 
gbflg(8) = 0 
gbflg(9) = 1 
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* note 

* 

Input angles are in array gbangles, input matrix is idmat, and 
rotation flags are in array gbflg 

*note 

Execute subfunction euler which performs angle rotation(s) 

* note 

Output matrix is gbmat 

* note 

* 

* 

In special cases, up to twelve knuckle angles may be used; 

the other nine knuckle angles are 

* 

* 

psij, thetaj, phij 

* 

* 

psip, thetap, phip 

* 

psiv, thetav, pliiv 


34 


Rotation of Angles from Gravity Axis to Model(Body) Axis 

Input Variables 

gbmat - 3x3 gravity to balance axis array (Euler input) 
gmangles - three-element array of angles (Euler input) 

gmflg - three-element array of rotation flags(Euler input); one flag for each angle 
where values are [- 1 = yaw rotation , 0 = pitch rotation, 1 = roll rotation ] 
phib - model(body) roll angle (deg) 
psib - model(body) yaw angle (deg) 
thetab - model(body) pitch angle (deg) 

Output Variables 

gmmat - 3x3 gravity to model(body) axis array (Euler output) 


* note gbmat and gmmat are double precision 

* note gmangles, gmflg, gbmat, and gmmat are not written to output datafile 


For a yaw-pitch-roll rotation, Let 

gmangles(l) = psib 
gmangles(2) = thetab 
gmangles(3) = phib 
gmflg(l) = -1 
gmflg(2) = 0 
gmflg(3) = 1 


* note Input angles are in array gmangles, input matrix is gbmat, and 

* rotation flags are in array gmflg 

*note Execute subfunction euler which performs angle rotation(s) 

* note Output matrix is gmmat 
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Rotation of Angles from Wind Axis to Gravity Axis 

Input Variables 

idmat - 3x3 identity array (Euler input) 
psiu - sideflow angle (deg) 
thetau - upflow angle (deg) 

wgangles - two-element array of angles (Euler input) 

wgflg - two-element array of rotation flags (Euler input); one flag for each angle 
where values are [0 = pitch rotation, - 1 = yaw rotation] 

Output Variables 

wgmat - 3x3 wind to gravity axis array (Euler output) 


* note idmat and wgmat are double precision 

* note idmat, wgflg, wgmat are not written to output datafile 

For a pitch-yaw rotation, Let 


wgangles(l) = thetau 
wgangles(2) = psiu 
wgflg(l) = 0 
wgflg(2) = -1 


* note 


Input angles are in array wgangles, input matrix is idmat, and 
rotation flags are in array wgflg 


*note 


Execute subfunction euler which performs angle rotation(s) 


* note 


Output matrix is wgmat 
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Rotation of Angles from Wind Axis to Model(Body) Axis 

Input Variables 

gmmat - 3x3 gravity to model(body) axis array 
wgmat - 3x3 wind to gravity axis array 

Output Variables 

wmmat - a 3x3 wind to model(body) axis array 

* note gmmat, wgmat, and wmmat are double precision 

* note gmmat, wgmat, and wmmat are not written to output datafile 

Let 


gmmat = 


and 


wgmat = 


then 


gmll gml2 gml3 

gm21 gm22 gm23 gravity to model(body) axis matrix 
gm31 gm32 gm33 


wgll wgl2 wgl3 

wg2 1 wg22 wg23 wind to gravity axis matrix 

wg31 wg32 wg33 


wmll wml2 wml3 

wmmat = wm21 wm22 wm23 wind to model(body) axis matrix 
wm31 wm32 wm33 


do k=l,3 
do j=l,3 

wm(j,k) = 0.0 
do i=l,3 

wm(j,k) = wm(j,k)+(gm(j,i)*wg(i,k)) 
enddo 
enddo 
enddo 
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Balance Axis Angles 

Input Variables 

gbmat - 3x3 gravity to balance axis array 
id - data identification flag 

Output Variables 

phi - balance roll angle[-360 to +360 degrees] (deg) 

phiO - balance roll angle[-360 to +360 degrees] at wind-off zero (deg) 

psi - balance yaw angle (deg) 

psiO - balance yaw angle at wind-off zero (deg) 

theta - balance pitch angle (deg) 

thetaO - balance pitch angle at wind-off zero (deg) 


* note rtod = 57.2957795131 (converts radians to degrees } 

* note gbmat is double precision 

* note gbmat, psi, and psiO are not written to output datafile 

* note datan2( value 1, value2) — datan( value 1 / value2) 

* where datan2 is double precision arc tangent of two arguments 


* gravity to balance axis angles 



gbl 1 

gbl2 

gbl 3 

gbmat = 

gb21 

gb22 

gb23 (gravity to balance axis matrix) 


gb31 

gb32 

gb33 


then 

if ( gbl2 = 0.0 & gbl 1 = 0.0 ) then 
psi = 0.0 
else 

psi = datan2(-gbl2, gbl 1) * rtod 
endif 

theta = dasin(-gbl3) * rtod 
if ( gb23 = 0.0 & gb33 = 0.0 ) then 
phi = 0.0 
else 

phi = datan2(-gb23, gb33) * rtod 
endif 
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note 


woz is wind-off zero flag number 


* 


if ( id = woz ) then 
psiO = psi 
thetaO = theta 
phiO = phi 
endif 


39 


Model(Body) Axis Angles 

Input Variables 


gmmat - 3x3 gravity to model(body) axis array 
wmmat - 3x3 wind to model(body)axis array 


Output Variables 

alfunc - uncorrected model(body) angle of attack (deg) 

beta - model(body) sideslip angle (deg) 

modrol - model(body) roll angle (deg) [-180 to +180 degrees; 

used for airplane models] 

pitchgml - model(body) pitch angle 
rollgml - model(body) roll angle 
yawgml - model(body) yaw angle 

* note rtod = 57.2957795131 (converts radians to degrees } 

* note gmmat and wmmat are double precision 

* note gmmat, wmmat, yawgml, pitchgml, and rollgml 

* are not written to output datafile 

* note datan2( value 1, value2) = datan( value 1 / value2) 

* where datan2 is double precision arc tangent of two arguments 


* gravity to model(body) axis angles 


Let 


gmmat = 


then 


gml 1 gml2 gml3 

gm21 gm22 gm23 gravity to model(body) axis matrix 
gm31 gm32 gm33 


if ( gml2 - 0.0 & gml 1 = 0.0 ) then 
yawgml = 0.0 
else 

yawgml = datan2(-gml2, gml 1) * rtod 
endif 

pitchgml = dasin(-gml3) * rtod 


if ( gm23 = 0.0 & gm33 = 0.0 ) then 
rollgml = 0.0 
else 

rollgml = datan2(-gm23, gm33) * rtod 
endif 
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* wind to model(body) axis angles 


Let 


wmmat 


l[k' I! 


wmll wml2 wml3 

wm21 wm22 wm23 wind to model(body) axis matrix 
wm31 wm32 wm33 


beta = -dasin(wm21) * rtod 

if ( wm3 1 = 0.0 & wml 1 = 0.0 ) then 
alfunc = 0.0 
else 

alfunc = datan2(wm31, wml 1) * rtod 
endif 

if ( wm23 = 0.0 & wm22 = 0.0 ) then 
modrol = 0.0 
else 

modrol = datan2(-wm23, wm22) * rtod 
endif 
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Model Height 

Input Variables 

elev - cart elevation plus cart offset (inches) 
pitchm - mast pitch angle (deg) 

rlength - distance from (1) the center of the vertical mast to the model height reference point or 

(2) the center of the pitch rotation pivot to the model height reference point; constant (inches) 
scode - sting code; constant 
thetas - strut pitch angle (deg) 

vlength - vertical distance from the pivot center on the vertical strut to the model 
reference point; constant (inches) 

Output Variables 

hgt - vertical distance from the test section floor to model height reference 
point (inches) 


* note dtor = 0.0174532925199 {converts degrees to radians} 


if ( scode = 0.0 ) then 

* no height calculation 

hgt = 0.0 

else if (( scode = 1.0 ) or ( scode = 2.0 )) then 

* height calculation for mast-mounted models or for alpha-beta sting 

hgt = 87.0+(elev-87.0)*cos(pitchm*dtor)+rlength*sin(pitchm*dtor) 
else if ( scode = 3.0 ) then 

* height calculation for vertical strut 

hgt = elev+rlength*sin(thetas*dtor)+vlength*cos(thetas*dtor) 
endif 
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Balance Components Rotated and Translated to Model Axis 

Input Variables 

afl.5 - axial force load corrected for method of attachment (lbs) 

bcmangles - three-element array of angles (Euler input) 

bcmflg - three-element array of rotation flags (Euler input); one flag for each 

angle where values are [ - 1 = yaw rotation, 0 = pitch rotation, 1 = roll rotation ] 
bcmin - six-element array (Euler input) 
bspanl - wing span; constant (ft) 
chord 1- wing aerodynamic chord; constant (ft) 
nfl.5 - normal force load corrected for method of attachment (lbs) 
phib - model(body) roll angle (deg) 

pm 1.5 - pitching moment load corrected for method of attachment (in-lbs) 
psib - model(body) yaw angle (deg) 

q - tunnel dynamic pressure corrected for compressibility; wall corrections (psf) 
rml.5 - rolling moment load corrected for method of attachment (in-lbs) 
sareal - wing area; constant (sqft) 

sfl.5 - side force load corrected for method of attachment (lbs) 
thetab - model(body) pitch angle (deg) 

xbarl - moment transfer distance(x direction) measured in the model force axis 
system from the moment center to the desired moment center(positive in 
the direction of positive model thrust(axial), side and normal force (ft) 
ybarl - moment transfer distance(y direction) measured in the model force axis 
system from the moment center to the desired moment center(positive in 
the direction of positive model thrust(axial), side and normal force (ft) 
yml.5 - yawing moment load corrected for method of attachment (in-lbs) 
zbarl - moment transfer distance(z direction) measured in the body force axis 

system from the moment center to the desired moment center(positive in 
the direction of positive model thrust(axial), side and normal force (ft) 

Output Variables 

bcmout - six-element array (Euler output) 

fa 1.1 - axial force rotated to the model axis (lbs) 

fa 1.2 - axial force rotated and translated to the model axis (lbs) 

fnl.l - normal force rotated to the model axis (lbs) 

fnl.2 - normal force rotated and translated to the model axis (lbs) 

fyl.l - side force rotated to the model axis (lbs) 

fyl.2 - side force rotated and translated to the model axis (lbs) 

mxl. 1 - rolling moment rotated to the model axis (in-lbs) 

mxl.2 - rolling moment rotated and translated to the model axis (in-lbs) 

my 1.1- pitching moment rotated to the model axis (in-lbs) 

my 1.2 - pitching moment rotated and translated to the model axis (in-lbs) 

mzl . 1 - yawing moment rotated to the model axis (in-lbs) 

mzl .2 yawing moment rotated and translated to the model axis (in-lbs) 
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* forces and moments rotated to the model axis 


For a yaw-pitch-roll rotation, Let 

bcmangles(l) = psib 
bcmangles(2) = thetab 
bcmangles(3) = phib 
bcmflg(l) = -1 
bcmflg(2) = 0 
bcmflg(3) = 1 

and 


bcmin(l) = afl.5 
bcmin(2) = sfl.5 
bcmin(3) = nfl.5 
bcmin(4) = -rml.5 
bcmin(5) = pml.5 
bcmin(6) = -yml.5 

* note execute subfunction euler which performs angle rotation(s) 

fa 1.1 = bcmout(l) 
fyl.l = bcmout(2) 
fnl.l = bcmout(3) 
mxl.l = -bcmout(4) 
my 1.1 = bcmout(5) 
mzl.l = -bcmout(6) 


* forces and moments rotated and translated to the model axis 


fal.2 = fal.l 
fyl.2 = fyl.l 
fill. 2 = fnl.l 

mxl.2 = mxl.l + ((fnl.l * ybarl) - (fyl.l * zbarl)) 
myl.2 = myl.l - ((fnl.l * xbarl) - (fal.l * zbarl)) 
mzl.2 = mzl.l - ((fyl.l * xbarl) - (fal.l * ybarl)) 
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Apply Estimates of Heyson’s Corrections - Walls Up 

Input Variables 

alfunc - uncorrected model(body) angle of attack (deg) 
dalup#.# - array of alpha correction coefficients; constants 
dolup# - array of drag over lift coefficients; constants 
dqup#.# - array of q correction coefficients; constants 
fa 1.2 - axial force rotated and translated to the model axis (lbs) 
fnl.2 - normal force rotated and translated to the model axis (lbs) 
nbal - balance number; constant 
pa - ambient pressure (psf) 

qu - uncorrected tunnel dynamic pressure with compressibility (psf) 
sareal - wing area; constant (sqft) 

velu - uncorrected test section free stream velocity (ft/sec) 

wcode - test section configuration code for classical wall corrections; constant 

Output Variables 

alpha - model(body) angle of attack with walls up or down corrections (deg) 

pstat - test section static pressure with walls up or down corrections (ft/sec) 

q - tunnel dynamic pressure corrected for compressibility; walls up or down corrections (psf) 

qcovrq - ratio of corrected dynamic pressure to uncorrected dynamic pressure 

vel - test section free stream velocity with walls up or down corrections (ft/sec) 

* note dtor = 0.0174532925199 { converts degrees to radians j 

* note clslu and cdslu are not written to output datafile 


* note 

* array dolup# is 


dolupl, dolup2, dolup3, dolup4, dolup5 


* array dalup#.# is 


dalup 1.1, dalup2.1, dalup3.1, 
dalupl.2, dalup2.2, dalup3.2, 
dalupl.3, dalup2.3, dalup3.3, 
dalupl.4, dalup2.4, dalup3.4, 
dalupl.5, dalup2.5, dalup3.5 


* 

* 

* 

* 

* 


* array dqup#.# is 


dqup 1.1, dqup2.1, dqup3.1, 
dqupl.2, dqup2.2, dqup3.2, 
dqupl.3, dqup2.3, dqup3.3, 
dqupl.4, dqup2.4, dqup3.4, 
dqupl.5, dqup2.5, dqup3.5 


* 

* 

* 

* 


* note 


execute walls up calculations if qu >0.1 nbal = 1 wcode = 1 
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* uncorrected stability axis lift coefficient 

clslu = ((fnl.2 * cos(alfunc * dtor)) - (fal.2 * sin(alfunc * dtor)))/ (qu * sareal) 

* uncorrected stability axis drag coefficient 

cdslu = ((fal.2 * cos(alfunc * dtor)) + (fnl.2 * sin(alfunc * dtor))) / (qu * sareal) 

* note execute subfunction walcor which computes corrected alpha and q 

* Tunnel parameters corrected for Q 

* note pstat calculation below is being reviewed 

vel = velu * sqrt(q / qu) 
pstat = pa - q 
qcovrq = q / qu 
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Apply Estimates of Heyson’s Corrections - Walls Down 

Input Variables 

alfunc - uncorrected model(body) angle of attack (deg) 
daldn#.# - array of alpha correction coefficients; constants 
doldn# - array of drag over lift coefficients; constants 
dqdn#.# - array of q correction coefficients; constants 
fa 1.2 - axial force rotated and translated to the model axis (lbs) 
fnl.2 - normal force rotated and translated to the model axis (lbs) 
nbal - balance number; constant 
pa - ambient pressure (psf) 

qu - uncorrected tunnel dynamic pressure with compressibility (psf) 
sareal - wing area; constant (sqft) 

velu - uncorrected test section free stream velocity (ft/sec) 

wcode - test section configuration code for classical wall corrections; constant 

Output Variables 

alpha - model(body) angle of attack with walls up or down corrections (deg) 

pstat - test section static pressure with walls up or down corrections (ft/sec) 

q - tunnel dynamic pressure corrected for compressibility; walls up or down corrections (psf) 

qcovrq - ratio of corrected dynamic pressure to uncorrected dynamic pressure 

vel - test section free stream velocity with walls up or down corrections (ft/sec) 

* note dtor = 0.0174532925199 (converts degrees to radians} 

* note clslu and cdslu are not written to output datafile 


* note 

* array doldn # is 


doldnl, doldn2, doldnS, doldn4, do!dn5 


* array daldn#.# is 


daldnl.l, daldn2.1, daldn3.1, 
daldnl.2, daldn2.2, daldn3.2, 
daldnl.3, daldn2.3, daldn3.3, 
daldnl.4, daldn2.4, daldn3.4, 
daldnl.5, daldn2.5, daldn3.5 


* 

* 

* 

* 

* 


* array dqdn#.# is 


dqdnl.l, dqdn2.1, dqdn3.1, 
dqdn 1.2, dqdn2.2, dqdn3.2, 
dqdnl.3, dqdn2.3, dqdn3.3, 
dqdn 1.4, dqdn2.4, dqdn3.4, 
clqdn 1.5, dqdn2.5, dqdn3.5 


* 

* 

* 

* 
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* note execute walls down calculations if qu > 0.1 nbal = 1 wcode = 2 

* uncorrected stability axis lift coefficient 

clslu = ((fnl.2 * cos(alfunc * dtor)) - (fal.2 * sin(alfunc * dtor)))/ (qu * sareal) 

* uncorrected stability axis drag coefficient 

cdslu = ((fal.2 * cos(alfunc * dtor)) + (fnl.2 * sin(alfunc * dtor))) / (qu * sareal) 

* note execute subfunction walcor which computes corrected alpha and q 

* Tunnel parameters corrected for Q 

* note pstat calculation below is being reviewed 

vel = velu * sqrt(q / qu) 
pstat = pa - q 
qcovrq = q / qu 
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Blockage and Jet Boundary Corrections 

Input Variables 

alfunc - uncorrected model(body) angle of attack (deg) 
beta - model(body) sideslip angle (deg) 

bodyblok - solid-blockage velocity effect for a body of revolution; constant 
bspanl - wing span; constant (ft) 

fa 1.2 - axial force rotated and translated to the model axis (lbs) 

fnl.2 - normal force rotated and translated to the model axis (lbs) 

jbcorr2 - jet boundary angle of attack factor; constant 

machu - uncorrected free stream mach number 

pstatu - uncorrected test section static pressure (psf) 

qu - uncorrected tunnel dynamic pressure with compressibility (psf) 

rhou - uncorrected air density 

rnftu - uncorrected test section reynolds number 

sareal - wing area; constant (sqft) 

tsarea - test section area; constant 

velu - uncorrected test section free stream velocity 

wcode - test section configuration code for classical wall corrections; constant 
wingblok - solid-blockage velocity effect for a wing; constant 

Output Variables 

alpha - model(body) angle of attack corrected for jet boundary (deg) 

delalp - jet boundary alpha correction 

mach - free stream mach number corrected for blockage 

pstat - test section static pressure; walls up or down or blockage corrections (psf) 
q - tunnel dynamic pressure corrected for compressibility; 

walls up or walls down or blockage corrections (psf) 
qcovrq - ratio of corrected dynamic pressure to uncorrected dynamic pressure 
rho - air density corrected for blockage (slugscf) 
rnft- test section reynolds number corrected for blockage (1/feet) 

vel - test section free stream velocity; walls up or walls down or blockage corrections (ft/sec) 


* note pie, csaf fdslu, clslu, cdlu, acd, machu2, and bcf are 

* not written to output datafile 

* note dtor = 0.0174532925199 {converts degrees to radians j 
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if (wcode = 3) then 
pie = 3.14159265 
csaf = sareal / (tsarea * 4.0) 

* uncorrected stability axis drag force 

fdslu = (fa 1.2 * cos(alfunc * dtor)) + (fnl.2 * sin(alfunc * dtor)) 

* uncorrected stability axis lift coefficient 

clslu = ((fnl.2 * cos(alfunc * dtor)) - (fal.2 * sin(alfunc * dtor))) / (qu * sareal) 

* uncorrected wind axis drag coefficient 

cdlu = ((fdslu * cos(beta * dtor)) - (fyl.2 * sin(beta * dtor))) / (qu * sareal) 

* Blockage calculations 

acd = cdlu - (clslu * clslu * sareal) / (pie * (bspanl / 12.0)**2) 
machu2 = machu * machu 

* jbcorr4(blockage correction) is wingblok + bodyblok 

bcf = (wingblok + bodyblok) / (1.0 - machu2)**1.5 + (1.0 + 0.4 * machu2) 

* acd * csaf / (1.0 - machu2) 

mach = machu * (1.0 + (1.0 + 0.2 * machu2) * bcf) 

pstat = pstatu * (1.0 - 1.4 * machu2 * bcf) 

q = qu * (1.0 + (2.0 - machu2) * bcf) 

qcovrq = q / qu 

rho = rhou * (1.0 - machu2 * bcf) 

rnft = rnftu * (1 .0 + (1.0 - 0.7 * machu2) * bcf) 

vel = velu * (1.0 + bcf) 

* Jet boundary angle of attack correction 

delalp = jbcorr2 * clslu 
alpha = alfunc + delalp 
endif 
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Heyson’s Boundary (Wall) Interference Functions 


Harry Heyson wrote several documents which described calculating boundary(wall) interference 
factors for a variety of configurations. In reference 2, Heyson discussed sixteen software functions. 
Curently, the 14- by 22-Foot Subsonic Tunnel’s software utilizes function 2. Any of the other functions 
may be implemented upon request. The function number and the type of interference for a specific model 
are as follows. 


Function 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 


Interference 

wind tunnel interference near a vanishingly small model 

average wind tunnel interference over a swept wing 

distribution od wind tunnel interference over the span of a swept wing 

average wind tunnel interference over a tail behind a swept wing 

average wind tunnel interference over a swept wing caused 

by the presence of lifting jets 

distribution of wind tunnel interference over the span of a 
swept wing caused by the presence of lifting jets 
average wind tunnel interference over a tail caused by the 
presence of lifting jets 

average wind tunnel interference over a single rotor 
distribution of wind tunnel interference over the lateral axis 
of a single rotor 

distribution of wind tunnel interference over the longitudinal 
axis of a single rotor 

average wind tunnel interference over a tail behind a single 
rotor 

average wind tunnel interference over tandem rotors 
average wind tunnel interference over unloaded rotor 
configuration 

average wind tunnel interference over a tail behind an 
unloaded rotor configuration 

average wind tunnel interference over side-by-side rotor 
configuration 
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16 


average wind tunnel interference over a tail behind a 
side-by-side rotor configuration 


Refer to references 2, 3, and 4 for details concerning Heyson’s boundary interference 
calculations. 
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Heyson’s Boundary (Wall) Interference 

Function 2 - Average Wind Tunnel Inteference Over a Swept Wing 

Input Variables 

alfunc - uncorrected model(body) angle of attack (deg) 
beta - model(body) sideslip angle (deg) 
bspanl - wing span; constant (ft) 

fa 1.2 - axial force rotated and translated to the model axis (lbs) 
fnl.2 - normal force rotated and translated to the model axis (lbs) 
gamma - ratio of tunnel width to height; constant 
hannon - used in etal calculation; constant 

hgt - height of the model height reference point above the test section floor (inches) 
lambda - angle of sweep angle; constant 
nbal - balance number; constant 

q - tunnel dynamic pressure with compressibility; walls up or walls down or blockage corrections (psf) 
rho - air density corrected for blockage (slugscf) 

rlength - distance from ( 1 ) the center of the vertical mast to the model height reference point 
or (2) the center of the pitch rotation pivot to the model height reference point; 
constant (inches) 

sareal - reference area; constant (sqft) 

semispan - semispan model flag; constant 

sigma - ratio of wing apn to the tunnel width; constant 

tsarea - test section area; constant (sqft) 

tsconf - test section configuration code for Hey son wall corrections; constant 
vlength - vertical distance from the pivot center on the vertical strut to the model reference point; 
constant (inches) 

wingar - wing aspect ratio; constant (ft) 
wingload - wind load configuration; constant 

zonecon - ratio of wind tunnel semi-height to height of model above tunnel floor; constant 
zoneflg - zetal calculation flag; constant 

Output Variables 

alfac - model(body) angle of attack with Heyson wall corrections (deg) 

cdc - drag coefficient with Heyson wall corrections 

clc - lift coefficient with Heyson wall corrections 

dragc - drag corrected with Heyson wall corrections 

liftc - lift corrected with Heyson wall corrections 

lode - lift over drag with Heyson wall corrections 

qc - tunnel dynamic pressure corrected for compressibility; blockage & Heyson wall corrections (psf) 
xeff - wake deflection angle (deg) 


* note dtor = 0.0174532925199 j converts degrees to radians} 

* note if hannon =0.0 rlength =0.0 vlength = 0.0 then etal =1.0 

* 14x22-F oot Subsonic Tunnel is 21.75 feet wide and 14.5 feet high 

* gamma = 21.75 / 14.5 = 1.5 

* sigma = bpsanl / 21.75 where bspanl is in feet 
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wingar = (bpsanl * bspanl) / sareal where bp sanl is in feet and 

sareal is in square feet 

tsarea = 21.75 * 14.5 = 315.375 

bhtw = half the tunnel width (inches); not written to output datafile 


bhtw = 10.875 * 12.0 
xeff = 0.0 

if (nbal =1.0 and tsconf > 0.0) then 

if ( zoneflg =1.0) then 
if ( hgt > 0.0 ) then 
zetal = zonecon / hgt 
else 

zetal = zonecon 
endif 
else 

zetal = zonecon 
endif 

xwid = (hannon+(rlength*cos(alpha*dtor))-(vlength*sin(alpha*dtor))) 
* sin(beta*dtor) 
etal = (bhtw-xwid)/bhtw 

if (semispan = 0.0) then 

* full wing span corrections 

nforce = fnl.2 
aforce = fa 1.2 
wrefarea = sareal 
tsecarea = tsarea 
wingspan = bspanl 

else if (semispan = 1.0) then 

* semi wing span corrections 

nforce = fnl.2 * 2.0 
aforce = fal.2 * 2.0 
wrefarea = sareal * 2.0 
tsecarea = tsarea * 2.0 
wingspan = bspanl * 2.0 

endif 
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* note execute subfunction hey son as follows 

* call hey son(tsconf, wingload, zetal, eta 1, gamma, sigma, lambda, 

* nforce, aforce, qoc, rhoc, alfunc, wrefarea, wingar, 

* tsecarea, wingspan, liftc, dragc, clc, cdc, qc, alfac, xeff) 

if (semispan = 1 .0) then 
liftc = liftc * 0.5 
dragc = dragc * 0.5 
endif 

lode = clc / cdc 
endif 


55 


Base(Cavity) Pressures 

Input Variables 

alpha - model angle of attack with jet boundary corrections (deg) 

arpba#.# - array of axial force base pressure areas; constant (sqft) 

arpbn#.# - array of normal force base pressure areas; constant (sqft) 

arpbp#.# - array of pitching moment base pressure moment arm*areas; constant(in-sqft) 

arpbr#.# - array of rolling moment base pressure moment arm*areas; constant (in-sqft) 

arpbs#.# - array of side force base pressure areas; constant (sqft) 

arpby#.# - array of yawing moment base pressure moment arm*areas; constant (in-sqft) 

bspanl - wing span; constant (ft) 

chord 1 - wing aerodynamic chord; constant (ft) 

phase# - array of base pressures (psf) 

pstat - tunnel static pressure with walls up/down or blockage corrections (psf) 

q - tunnel dynamic pressure with compressibility; walls up or walls down or blockage corrections (psf) 
sareal - wing area; constant (sqft) 

Output Variables 

cabase 1 - axial force base pressure tare coefficient calculated using axial force 
base pressure tare applied in final axial force calculation 
cdbasel - drag base pressure tare coefficient 
cnbasel - normal force base pressure tare coefficient 
cpbase# - array of base pressure coefficients 
cpmbasel - pitching moment base pressure tare coefficient 
crmbasel - rolling moment base pressure tare coefficient 
cybasel - side force base pressure tare coefficient calculated using side force 
base pressure tare applied in final side force calculation 
cymbasel - yawing moment base pressure tare coefficient 
dpbase# - array of differential base pressures (psf) 
fabasel - axial force base pressure tare (lbs) 
fnbasel - normal force base pressure tare (lbs) 
fybasel - side force base pressure tare (lbs) 
pmbase 1 - pitching moment base pressure tare (ft-lbs) 
rmbase 1 - rolling moment base pressure tare (ft-lbs) 
ymbase 1 - yawing moment base pressure tare (ft-lbs) 

* note maximum number ofbase( cavity) pressures per balance is 25 

* cavity pressures are calculated using base pressure names & equations 

* base(cavity) pressure areas (constants with default values of 0.0) 

arpba 1.1 (axial force; base pressure 1, balance 1) 

arpba2.1 (axial force; base pressure 2, balance 1) 

arpba3.1 (axial force; base pressure 3, balance 1) 


arpba23.1 


(axial force; base pressure 23, balance 1) 
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arpba24. 1 (axial force; base pressure 24, balance 1) 

arpba25.1 (axial force; base pressure 25, balance 1) 


arpbyl.l 
arpby2. 1 
arpby3.1 


(yawing moment; base pressure 1, balance 1) 
(yawing moment; base pressure 2, balance 1) 
(yawing moment; base pressure 3, balance 1) 


arpby23.1 
arpby24. 1 
arpby25.1 


(yawing moment; base pressure 23, balance 1) 
(yawing moment; base pressure 24, balance 1) 
(yawing moment; base pressure 25, balance 1) 


* differential base(cavity) pressures 


dpbasel = pbasel-pstat 
dpbase2 = pbase2-pstat 
dpbase3 = pbase3-pstat 


( base pressure 1 ) 
( base pressure 2 ) 
( base pressure 3 ) 


dpbase23 = pbase23-pstat (base pressure 23) 

dpbase24 = pbase24-pstat ( base pressure 24 ) 

dpbase25 = pbase25-pstat (base pressure 25) 

* base(cavity) pressure coefficients 


cpbasel = dpbasel /q 
cpbase2 = dpbase2/q 
cpbase3 = dpbase3/q 


( base pressure 1 ) 
( base pressure 2 ) 
( base pressure 3 ) 


cpbase23 = dpbase23/q ( base pressure23 ) 

cpbase24 = dpbase24/q ( base pressure 24 ) 

cpbase25 = dpbase25/q ( base pressure 25 ) 

* base(cavity) pressure tares 

fabasel= -( (dpbasel*arpbal.l) + (dpbase2*arpba2.1) + (dpbase3*arpba3.1) + .... 

+ (dpbase23*arpba23.1) + (dpbase24*arpba24.1)+ (dpbase25*arpba25.1) ) 

fybasel = (dpbasel*arpbsl.l) + (dpbase2*arpbs2.1)+ (dpbase3*arpbs3.1) + .... 

+ dpbase23*arpbs23.1) + (dpbase24*arpbs24.1)+ (dpbase25*arpbs25.1) 

fnbasel = (dpbasel*arpbnl.l) + (dpbase2*arpbn2.1)+ (dpbase3*arpbn3.1) + .... 

+ dpbase23*arpbn23.1) + (dpbase24*aipbn24. 1)+ (dpbase25*arpbn25.1) 
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rmbasel = (dpbasel*arpbrl.l) + (dpbase2*arpbr2.1)+ (dpbase3*arpbr3.1) + .... 

+ dpbase23*arpbr23.1) + (dpbase24*arpbr24. 1)+ (dpbase25*arpbr25.1) 

pmbasel = (dpbasel*arpbpl.l) + (dpbase2*arpbp2.1)+ (dpbase3*arpbp3.1) + .... 

+ dpbase23*arpbp23.1) + (dpbase24*arpbp24. 1)+ (dpbase25*arpbp25.1) 

ymbasel = (dpbasel*arpbyl.l) + (dpbase2*arpby2.1)+ (dpbase3*arpby3.1) + .... 

+ dpbase23*arpby23.1) + (dpbase24*arpby24. 1)+ (dpbase25*arpby25.1) 


* note dtor = 0.0174532925199 {converts degrees to radians} 

* base(cavity) pressure tare coefficients 

cabasel = fabasel/(q*sareal) 
cybasel = fybasel/(q*sareal) 
cnbasel = fnbasel/(q*sareal) 
crmbasel = rmbasel/(q*sareal*bspanl) 
cpmbasel = pmbasel/(q*sareal*chordl) 
cymbasel = ymbasel/(q*sareal*bspanl) 

* base(cavity) pressure tare drag coefficient 

cdbasel = (cabasel *cos(alpha*dtor))+(cnbasel*sin(alpha*dtor)) 


58 


Model Axis Components and Coefficients 

Input Variables 

bspanl - wing span; constant (ft) 

chord 1- wing aerodynamic chord; constant (ft) 

fabasel - axial force base pressure tare (lbs) 

fnbasel - normal force base pressure tare (lbs) 

fybasel - side force base pressure tare (lbs) 

pmbase 1 - pitching moment base pressure tare (ft-lbs) 

q - tunnel dynamic pressure corrected for compressibility; walls up or walls down or 
blockage corrections (psf) 
rmbase 1 - rolling moment base pressure tare (ft-lbs) 
sareal - wing area; constant (sqft) 
ymbase 1 - yawing moment base pressure tare (ft-lbs) 

Output Variables 

cal - model axis axial force coefficient corrected for base pressure tare 

cmxl - model axis rolling moment coefficient corrected for base pressure tare 

cmy 1 - model axis pitching moment coefficient corrected for base pressure tare 

cmzl - model axis yawing moment coefficient corrected for base pressure tare 

cnl - model axis normal force coefficient corrected for base pressure tare 

cy 1 - model axis side force coefficient corrected for base pressure tare 

fal - model(body) axis axial force corrected for base pressure tare (lbs) 

fa 1.2 - axial force rotated and translated to the model axis (lbs) 

fnl - model(body) axis normal force corrected for base pressure tare (lbs) 

fnl.2 - normal force rotated and translated to the model axis (lbs) 

fyl - model(body) axis side force corrected for base pressure tare (lbs) 

fyl.2 - side force rotated and translated to the model axis (lbs) 

mxl - model(body) axis rolling moment corrected for base pressure tare (in-lbs) 

mxl.2 - rolling moment rotated and translated to the model axis (in-lbs) 

my 1 - model(body) axis pitching moment corrected for base pressure tare (in-lbs) 

my 1.2- pitching moment rotated and translated to the model axis (in-lbs) 

mzl - model(body) axis yawing moment corrected for delcm and base pressure tare (in-lbs) 

mzl.2 - yawing moment rotated and translated to the model axis (in-lbs) 


* balance components 

fal = fal. 2 
fyl = fyl.2 
fnl = fnl.2 
mxl = mxl.2 
myl = myl.2 
mzl = mzl.2 
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* balance components & coefficients corrected for base pressure tares 


fal = fal - fabasel 
fyl = fyl - fybasel 
fnl = fnl - fnbasel 
mxl = mxl - rmbasel 
myl = myl - pmbasel 
mzl = mzl - ymbasel 

cal = fal / (q * sareal) 

cyl = fyl / (q * sareal) 

cnl = fnl / (q * sareal) 

cmxl = mxl / (q * sareal * bspanl) 

cmyl = myl / (q * sareal * chordl) 

cmzl = mzl / (q * sareal * bspanl) 
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Stability Axis Components and Coefficients 

Input Variables 

alpha - model angle of attack with jet boundary corrections (deg) 

bspanl - wing span; constant (ft) 

chord 1 - wing aerodynamic chord; constant (ft) 

fal - model(body) axis axial force corrected for base pressure tare (lbs) 
fnl - model(body) axis normal force corrected for base pressure tare (lbs) 
fyl - model(body) axis side force corrected for base pressure tare (lbs) 
mxl - model(body) axis rolling moment corrected for base pressure tare (in-lbs) 
my 1 - model(body) axis pitching moment corrected for base pressure tare (in-lbs) 
mzl - model(body) axis yawing moment corrected for base pressure tare (in-lbs) 
q - tunnel dynamic pressure corrected for compressibility; walls up or walls down or 
blockage corrections (psf) 
sareal - wing area; constant (sqft) 

Output Variables 

cds 1 - stability axis drag coefficient corrected for base pressure tare 

clsqrl - lift coefficient squared corrected for base pressure tare 

els 1 - stability axis lift coefficient corrected for base pressure tare 

cmxs 1 - stability axis rolling moment coefficient corrected for base pressure tare 

cmys 1 - stability axis pitching moment coefficient corrected for base pressure tare 

cmzsl - stability axis yawing moment coefficient corrected for base pressure tare 

cys 1 - model axis side force coefficient corrected for base pressure tare 

dell - jet boundary drag correction 

delp - jet boundary pitch correction 

fds 1 - stability axis drag corrected for base pressure tare (lbs) 

fls 1 — stability axis lift corrected for base pressure tare (lbs) 

fys 1 - stability axis side force corrected for base pressure tare (lbs) 

ls/ds 1 - lift-over-drag ratio corrected for base pressure tare 

mxs 1 - stability axis rolling moment corrected for base pressure tare (in-lbs) 

mys 1 - stability axis pitching moment corrected for base pressure tare (in-lbs) 

mzs 1 - stability axis yawing moment corrected for base pressure tare (in-lbs) 


* note dtor = 0.0174532925199 {converts degrees to radians} 

* balance components & coefficients corrected for base pressure tares 

fdsl = (fal * cos(alpha*dtor)) + (fnl * sin(alpha*dtor)) 
fysl = fyl 

fls 1 = (fnl * cos(alpha*dtor)) - (fal * sin(alpha*dtor)) 
mxsl = (mxl * cos(alpha*dtor)) + (mzl * sin(alpha*dtor)) 
mysl = myl 

mzsl = (mzl * cos(alpha*dtor)) - (mxl * sin(alpha*dtor)) 
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* jet boundary drag force correction 

dell = (jbcorrl * flsl * flsl) / (q * sareal) 
fdsl = fdsl + dell 

* jet boundary pitching moment correction 

delp = jbcorr3 * flsl * chord 1 
mys 1 = mys 1 - delp 

cdsl = fdsl / (q * sareal) 

cysl = fysl / (q * sareal) 

clsl = flsl / (q * sareal) 

cmxsl = mxsl / (q * sareal * bspanl) 

cmysl = mysl / (q * sareal * chord 1) 

cmzsl = mzsl / (q * sareal * bspanl) 

ls/dsl = clsl / cdsl 
clsqrl = clsl * clsl 
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Wind Axis Components and Coefficients 

Input Variables 

beta - model sideslip angle (deg) 
bspanl - wing span; constant (ft) 
chord 1 - wing aerodynamic chord; constant (ft) 
fds 1 - stability axis drag corrected for base pressure tare (lbs) 
fls 1 — stability axis lift corrected for base pressure tare (lbs) 
fys 1 - stability axis side force corrected for base pressure tare (lbs) 
mxs 1 - stability axis rolling moment corrected for base pressure tare (in-lbs) 
mys 1 - stability axis pitching moment corrected for base pressure tare (in-lbs) 
mzs 1 - stability axis yawing moment corrected for base pressure tare (in-lbs) 
q - tunnel dynamic pressure corrected for compressibility; walls up or walls down or 
blockage corrections (psf) 
sareal - wing area; constant (sqft) 

Output Variables 

cc 1 - wind axis crosswind coefficient corrected for base pressure tare 

cdl - wind axis drag coefficient corrected for base pressure tare 

ell - wind axis lift coefficient corrected for base pressure tare 

cmxwl - wind axis rolling moment coefficient corrected for base pressure tare 

cmywl - wind axis pitching moment coefficient corrected for base pressure tare 

cmzw 1 - wind axis yawing moment coefficient corrected for base pressure tare 

fc 1 - wind axis crosswind corrected for base pressure tare (lbs) 

fdl - wind axis drag corrected for base pressure tare (lbs) 

fll — wind axis lift corrected for base pressure tare (lbs) 

ldl - lift-over-drag ratio corrected for base pressure tare 

mxwl - wind axis rolling moment corrected for base pressure tare (in-lbs) 

mywl - wind axis pitching moment corrected for base pressure tare (in-lbs) 

mzw 1 - wind axis yawing moment corrected for base pressure tare (in-lbs) 


* note dtor = 0.0174532925199 {converts degrees to radians} 


* balance components & coefficients corrected for base pressure tares 

fdl = (fdsl * cos(beta*dtor)) - (fysl * sin(beta*dtor)) 
fcl = (fysl * cos(beta*dtor)) + (fdsl * sin(beta*dtor)) 
fll = flsl 

mxwl = (mxsl * cos(beta*dtor)) + (mysl * sin(beta*dtor)) 
mywl = (mysl * cos(beta*dtor)) - (mxsl * sin(beta*dtor)) 
mzw 1 = mzs 1 
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cdl = fdl / (q * sareal) 

ccl = fcl / (q * sareal) 

ell = fll / (q * sareal) 

cmxwl = mxwl / (q * sareal * bspanl) 

cmywl = mywl / (q * sareal * chord 1) 

cmzwl = mzwl / (q * sareal * bspanl) 

ldl = ell / cdl 
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Reference Axis Components and Coefficients 

Input Variables 

bspanl - wing span; constant (ft) 

chord 1- wing aerodynamic chord; constant (ft) 

fal - model(body) axis axial force corrected for base pressure tare (lbs) 

fnl - model(body) axis normal force corrected for base pressure tare (lbs) 

fyl - model(body) axis side force corrected for base pressure tare (lbs) 

mxl - model(body) axis rolling moment corrected for base pressure tare (in-lbs) 

my 1 - model(body) axis pitching moment corrected for base pressure tare (in-lbs) 

mzl - model(body) axis yawing moment corrected for base pressure tare (in-lbs) 

phirl - reference axis roll angle (deg) 

psirl - reference axis yaw angle (deg) 

q - tunnel dynamic pressure corrected for compressibility; walls up or walls down or 
blockage corrections (psf) 

refangles - three-element array of angles (Euler input) 
refflg - three-element array of rotation flags (Euler input); 

one flag for each angle where values are [ - 1 = yaw rotation, 0 = pitch rotation, 1 = roll rotation ] 
refin - six-element array (Euler input) 
sareal - wing area; constant (sqft) 
thetarl - reference axis pitch angle (deg) 

xref 1 - moment transfer distance(x direction) measured in the model force axis system from the 
moment center to the desired moment center(positive in the direction of positive model 
thrust( axial), side and normal force (ft) 

yref 1 - moment transfer distance(y direction) measured in the model force axis system from the 
moment center to the desired moment center(positive in the direction of positive model 
thrust( axial), side and normal force (ft) 

zrefl - moment transfer distance(z direction) measured in the body force axis system from the 
moment center to the desired moment center(positive in the direction of positive model 
thrust( axial), side and normal force (ft) 

Output Variables 

faref 1 - model(body) axis axial force rotated and translated to an arbitrary reference axis (lbs) 

farefrl - model(body) axis axial force rotated to an arbitrary reference axis (lbs) 

fnrefl - model (body) axis normal force rotated and translated to an arbitrary reference axis (lbs) 

fnrefrl - model (body) axis normal force rotated to an arbitrary reference axis (lbs) 

fyrefl - model (body) axis side force rotated and translated to an arbitrary reference axis (lbs) 

fyrefrl - model (body) axis side force rotated to an arbitrary reference axis (lbs) 

mxref 1 - model (body) axis rolling moment rotated and translated to an arbitrary reference axis (in-lbs) 
mxrefrl - model (body) axis rolling moment rotated to an arbitrary reference axis (in-lbs) 
myref 1 - model (body) axis pitching moment rotated and translated to an arbitrary reference 
axis (in-lbs) 

myrefrl - model (body) axis pitching moment rotated to an arbitrary reference axis (in-lbs) 

mzrefl - model (body) axis yawing moment rotated and translated to an arbitrary reference axis (in-lbs) 

mzrefrl - model (body) axis yawing moment rotated to an arbitrary reference axis (in-lbs) 

refout - six-element array (Euler output) 
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* forces and moments rotated to the model axis 


For a yaw-pitch-roll rotation, Let 

refangles(l) = psirl 
refangles(2) = thetarl 
refangles(3) = phirl 
refflg(l) = -l 
refflg(2) = 0 
refflg(3) = 1 


and 


refin(l) = fal 
refin(2) = fy 1 
refin(3) = fnl 
refin(4) = -mxl 
refin(5) - myl 
refin(6) = -mzl 

* note execute subfunction euler which performs angle rotation(s) 

* model(body) axis components rotated to an arbitrary reference axis 

farefrl =refout(l) 
fyrefrl = refout(2) 
fnrefrl = refout(3) 
mxrefrl = -refout(4) 
myrefrl = refout(5) 
mzrefrl = -refout(6) 

* model(body) axis components and coefficients rotated and translated to 

* an arbitrary reference axis 

farefl = farefrl 
fyrefl = fyrefrl 
fnrefl = fnrefrl 

mxrefl = mxrefrl + (fnrefrl * yrefl) - (fyrefrl * zrefl) 
myrefl = myrefrl - (fnrefrl * xrefl) - (farefrl * zrefl) 
mzrefl = mzrefrl - (fyrefrl * xrefl) - (farefrl * yrefl) 

carefl = farefl / (q * sareal) 
cyrefl = fyrefl / (q * sareal) 
cnrefl = fnrefl / (q * sareal) 
cmxrefl = mxrefl / (q * sareal * bspanl) 
cmyrefl = myrefl / (q * sareal * chordl) 
cmzrefl = mzrefl / (q * sareal * bspanl) 
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Model and Wall Pressures 


Input Variables 

dpat - difference between the settling chamber total pressure and the ambient pressure “ptot - pa” (psf) 
dpinf - dynamic pressure uncorrected for compressibility (psf) 
dprt - WICS set pressure (psf) 
p# - array of esp pressures (psf) 

q- tunnel dynamic pressure corrected for compressibility; wall corrections (psf) 
qu - uncorrected tunnel dynamic pressure with compressibility (psf) 
wsp# - array of WICS esp pressures (psf) 

Output Variables 

cpp# - array of pressure coefficients 
wpp# - array of WICS pressure coefficients 

* note maximum number of esp modules = 16 

* maximum number of pressures per esp module - 64 

* note qcalc is not written to output datafile 

* pressure coefficients 

if (q.eq.0.0) then 
qcalc = 1.0 
else 

qcalc = q 
endif 

cpplOl = (pl01-dpat+dpinf)/qcalc 
cpp 102 = (pl02-dpat+dpinf)/qcalc 
cpp 103 = (pl03-dpat+dpinf)/qcalc 


cpp 162 = (pl62-dpat+dpinf)/qcalc 
cpp 163 = (pl63-dpat+dpinf)/qcalc 
cpp 164 = (pl64-dpat+dpinf)/qcalc 


(module 1, port 1) 
( module 1, port 2 ) 
( module 1, port 3 ) 


(module 1, port 62) 
( module 1, port 63 ) 
(module 1, port 64) 


cpp 1601 = (p 1 60 1 -dpat+dpinf)/qcalc ( module 16, port 1 ) 
cpp 1602 = (pl602-dpat+dpinf)/qcalc ( module 16, port 2 ) 
cpp 1603 = (pl603-dpat+dpinf)/qcalc (module 16, port 3) 


cpp 1 662 = (p 1 662-dpat+dpinf)/qcalc ( module 16, port 62 ) 
cpp 1 663 = (p 1 663-dpat+dpinf)/qcalc ( module 16, port 63 ) 
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cpp 1 664 = (p 1 664-dpat+dpinf)/qcalc ( module 16, port 64 ) 


* Wall Interference Correction System pressures 

if (qu.eq.0.0) then 
qcalc =1.0 
else 

qcalc = qu 
endif 


South Wall Pressures 

wpp5001 = (wsp5001+dpinf-dprt)/qcalc 
wpp5002 = (wsp5002+dpinf-dprt)/qcalc 
wpp5003 = (wsp5003+dpinf-dprt)/qcalc 


wpp5062 = (wsp5062+dpinf-dprt)/qcalc 
wpp5063 = (wsp5063+dpinf-dprt)/qcalc 
wpp5064 = (wsp5064+dpinf-dprt)/qcalc 


wpp5 101 = (wsp5101 +dpinf-dprt)/qc ale 
wpp5102 = (wsp5102+dpinf-dprt)/qcalc 
wpp5103 = (wsp5103+dpinf-dprt)/qcalc 


wpp5162 = (wsp5162+dpinf-dprt)/qcalc 
wpp5163 = (wsp5163+dpinf-dprt)/qcalc 
wpp5164 = (wsp5164+dpinf-dprt)/qcalc 


Ceiling Pressures 

wpp5201 = (wsp5201+dpinf-dprt)/qcalc 
wpp5202 = (wsp5202+dpinf-dprt)/qcalc 
wpp5203 = (wsp5203+dpinf-dprt)/qcalc 


wpp5262 = (wsp5262+dpinf-dprt)/qcalc 
wpp5263 = (wsp5263+dpinf-dprt)/qcalc 
wpp5264 = (wsp5264+dpinf-dprt)/qcalc 


wpp5301 = (wsp5301+dpinf-dprt)/qcalc 
wpp5302 = (wsp5302+dpinf-dprt)/qcalc 
wpp5303 = (wsp5303+dpinf-dprt)/qcalc 


( module 50, port 1 ) 
( module 50, port 2 ) 
(module 50, port 3) 


( module 50, port 62 ) 
( module 50, port 63) 
(module 50, port 64) 

( module 51, port 1 ) 
(module 51, port 2) 
(module 51, port 3) 

(module 51, port 62) 
(module 51, port 63) 
(module 51, port 64) 


( module 52, port 1 ) 
( module 52, port 2 ) 
( module 52, port 3 ) 


( module 52, port 62 ) 
( module 52, port 63 ) 
( module 52, port 64 ) 
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( module 53, port 1 ) 
( module 53, port 2 ) 
(module 53, port 3) 


wpp5362 = (wsp5362+dpinf-dprt)/qcalc 
wpp5363 = (wsp5363+dpinf-dprt)/qcalc 
wpp5364 = (wsp5364+dpinf-dprt)/qcalc 


North Wall Pressures 

wpp5401 = (wsp5401+dpinf-dprt)/qcalc 
wpp5402 = (wsp5402+dpinf-dprt)/qcalc 
wpp5403 = (wsp5403+dpinf-dprt)/qcalc 


wpp5462 = (wsp5462+dpinf-dprt)/qcalc 
wpp5463 = (wsp5463+dpinf-dprt)/qcalc 
wpp5464 = (wsp5464+dpinf-dprt)/qcalc 


wpp5501 = (wsp5501+dpinf-dprt)/qcalc 
wpp5502 = (wsp5502+dpinf-dprt)/qcalc 
wpp5503 = (wsp5503+dpinf-dprt)/qcalc 


wpp5562 = (wsp5562+dpinf-dprt)/qcalc 
wpp5563 = (wsp5563+dpinf-dprt)/qcalc 
wpp5564 = (wsp5564+dpinf-dprt)/qcalc 


South Wall Centerline Pressures 

wpp5701 = (wsp5701+dpinf-dprt)/qcalc 
wpp5702 = (wsp5702+dpinf-dprt)/qcalc 
wpp5703 = (wsp5703+dpinf-dprt)/qcalc 


wpp5730 = (wsp5730+dpinf-dprt)/qcalc 
wpp5731 = (wsp5731+dpinf-dprt)/qcalc 
wpp5732 = (wsp5732+dpinf-dprt)/qcalc 


Ceiling Centerline Pressures 

wpp5801 = (wsp5801+dpinf-dprt)/qcalc 
wpp5802 = (wsp5802+dpinf-dprt)/qcalc 
wpp5803 = (wsp5803+dpinf-dprt)/qcalc 
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( module 53, port 62 ) 
( module 53, port 63 ) 
( module 53, port 64 ) 


( module 54, port 1 ) 
( module 54, port 2 ) 
( module 54, port 3) 


( module 54, port 62 ) 
( module 54, port 63) 
( module 54, port 64 ) 


( module 55, port 1 ) 
( module 55, port 2 ) 
( module 55, port 3 ) 


( module 55, port 62 ) 
( module 55, port 63 ) 
( module 55, port 64 ) 


(module 57, port 1 ) 
(module 57, port 2) 
(module 57, port 3) 


(module 57, port 30) 
(module 57, port 31) 
( module 57, port 32) 


( module 58, port 1 ) 
( module 58, port 2 ) 
(module 58, port 3) 


wpp5830 = (wsp5830+dpinf-dprt)/qcalc (module 58, port 30) 

wpp583 1 = (wsp583 l+dpinf-dprt)/qcalc f module 58, port 31 ) 

wpp5832 = (wsp5832+dpinf-dprt)/qcalc (module 58, port 32) 


North Wall Centerline Pressures 

wpp5901 = (wsp5901+dpinf-dprt)/qcalc (module 59, port 1) 

wpp5902 = (wsp5902+dpinf-dprt)/qcalc ( module 59, port 2 ) 

wpp5903 = (wsp5903+dpinf-dprt)/qcalc (module 59, port 3) 


wpp5930 = (wsp5930+dpinf-dprt)/qcalc (module 59, port 30) 

wpp593 1 = (wsp593 l+dpinf-dprt)/qcalc f module 59, port 31 ) 

wpp5932 = (wsp5932+dpinf-dprt)/qcalc (module 59, port 32) 
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